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EXECUTIVE SUMMAU'S 


This PinaX Report (in Tive volumes) was prepared hy 
Ocean Data Systems, Inc, tor the Jet Propulsion Laboratory 
under Contract Ho, 954668 (as amended) in support of the 
SEAShT Program, 

The objective was to develop atmospheric analysis and 
prediction models of varying (grid) resolution, and to test 
the models vising real observational data for the purpose of 
assessing the impact of grid resolution on short-range 
numerical weather prediction. Secondly, tlve discretionary 
model procedures had to be examined so that the computational 
viability of SEASAT data might be enhanced during the conduct 
of (future) sensitivity tests. 

ODSI has: (1) completed and tested an entire sequence of 
objective analysis models; (2) developed and tested four ver- 
sions of an atmospheric prediction model on grids of differing 
horizontal and/or vortical resolution; (3) prepared a real 
data base consisting of two four-day meteorological scenarios 
(using the analysis models) ; (4) exercised all analysis and 
prediction models (on the CYBER 175 at PNWC and the CDC 7600 
at NASA Ames) ; and (5) conducted sensitivity studies with the 
forecast medal to order the impacts of individual procedures 
and effects. 
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A, ' Genera.1 Descriptiion 

ThesD hemispheric models were designed for use on either 
of two horizontal grids: (l) a 63 x S3 coarse-inesh grid 
(CMG) ; and (2) a 187 x 187 fine-mesh grid (3?MG) , On the 
CMGf the true-earth mesh distance varies from 204 kilometers 
at the equator to 408 kilometers at the pole. On the ITMGf 
it varies from G8 to 136 kilometers, 5?he scales of meteoro- 
logical interest can be handled very well on the PMG, 

The objective analysis models are based on a Pattern 
Conserving Technique (PGT) which accommodates observations 
while preserving differential properties of the first-guess 
field (under a system of weights). In tlie multi-level upper- 
air analysis, the 1000- to 100-MB region was resolved with 
ten layers for better representation of stable layers in the 
lower troposphere. Stratospheric ejctensions (from 100 to 
10 MBS) to the mass structure were obtained using regression 
equations. Although analyses were performed on the FMG, the 
forecast models were initialised using interpolated CMG values. 

In the analysis effort, ODSIi (1) examined procedures for 
allowing data to influence the analysis f (2) examined the 
effects of varying the x^eights in the analysis procedure ? 

(3) tested and implemented procedures for solving the minimisation 
equation in an optimal x/ayr (4) described the impact of grid 
resolxxticn on analysis? and (S) devised and implemented 
numerous practical solutions to analysis problems, generally. 
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Wio IiQinlaphcrlG pjrotUction modols contain suitable 
numo3.'ios and physics for genei’ating operationally-compGtitivG 
forecasts, Bach model was cKercised in either a 5--Xayor or 
lQ--layoi: mode on either the CHG or PMG, echo contribution 
of each process f proGodure and. paratnoter to a forecast was 
ordered and oxaminod. Clearly , the primary consideration 
in prediction is to pose the initial conditions properly 
(through adequate data, resolution and initialisation methods) . 
Complex forecast models, however, employ "computational 
devices'' (such as smoothers, tendency trunoators, mountain 
ahapiiUTr diffusion operators, temporal filters) for both 
stabi3.ity and cosmetic purposes. The use/misuso of such 
devices can mask the effects of modeled processes and effects. 

The effects of resolution on prediction modeling are 
presented in terms of: (1) forecast charts; (2) forecast 
change/o: ror charts? (3) forecast difference charts (10-layer 
minus 5-layor or 187 s 187 forecast minus 63 x 63 forecast) ; 
and ('IN statistical/graphical information describing model 
performance and model behavior (energetics) , 

3, Sximmary of Findings 

Increases in horiaontal resolution arc computationally 
expensive, but they load te significant improvements In model 
performance and boJiavior, For a given number of vertical 
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layers ; lilio three-foM .tnerease in horinonlsal resolution 
represents a 27K problem* Numorioally/ the benefits are; 

(1) wore realistic phase speet^s of fronts/ troughs? (2) better 
representation of graslionts? (3) the possibility of treating 
attditional smaller scales? and (4) the redviction of harmful 
effects of computational devices* Xn terms of the physics ^ 
the benefits are; (1) the scale of dissipation becomes further 
removed from those scales producing eddy hinetic energy? (2) 
the effects of parameteriaed physical processes are. better 
retained in small scales r instead of being smeared upscale? 

(3) a 100?; incj^ease in modeled larocipitation; and (4) bettor 
treatment of energy conversions in baroelinic scales (cyelO” 
genesis) ♦ 

Ir.e-r eases in vertical resolution tend to produce smaller 
effects, 'the tcn<-*layer models pjrcdicted 10"20'i less precipi- 
tation than the five-layer models (on both the PMG and CMG) . 
There was no significant improvement indicated with respect 
to prediction of pressure system intensity. Increased ver- 
tical resolution permits improved specification of high- 
altitude conditions# but the impact on sea-level pressure 
forecasts is minimal. 

Analysis models must be capable of distributing obsor%»ed 
informat Ion into scales which are computa t i onally viable in 
the prediction models. More study of data influence functions 
and dynamic coupling of atmospheric variables is needed. 
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In the prediction models r the need for dynamic initiali- 
sation is readily apparent, such procedures may be compu- 
tationally expensive, but they are needed to: (1) minimise 
the harmful effects of initialisation shock? (2) improve 
model energetics? axid (3) Mtabilise the model context to 
prevejxt "wash out" of information in the smaller scales 
(such as SEAShT data) . 

Finally, the effects of discretionary procedures in 
prediction models must be understood, ordered and exploited. 
Physical effects are often small compared to the effects of 
" cosmetic operators " . 
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INTRODUCTION 


This is the Pinal Report of Ocean Data Systems, Inc. 
(ODSI) for the Jet Propulsion Laboratory (JPL) under Contract 
No. 954668 in support of the SEASAT Program of the National 
Aeronautics and Space Administration (NASA) . The Final 
Report is in five volumes, and describes: (1) atmospheric 
analysis and prediction models for grids of varying reso- 
lution; (2) model test results when executed using real data 
scenarios; and (3) , sensitivity studies with a coarse-raesh 
prediction model. 

Volxime I is the Summary of Findings for all aspects of 
the work defined in the original Statement of Work; i.e., 
the development and testing of analysis models and the first 
three (of four) versions of a prediction model. Volumes II 
and III contain rather detailed descriptions of the analysis 
and prediction models, respectively. Volume IV summarizes 
the results of sensitivity studies conducted with a coarse- 
mesh prediction model. Volume V contains an update of 
portions of Volume I to reflect the tests conducted with a 
fourth version of the prediction model (the 187 x 187 x 10 
layer model) . 



The NASA-JPL-ODSI atmospheric modeling effort derives 
from the need to develop and test an operationally-oriented 
model capability for assessing the impact of SEASAT data on 
objective analysis and short-range (a few days) numerical 
weather prediction. As a practical matter, it means that a 
proper model context has to be assembled which facilitates 
the necessary sensitivity studies. Problem areas {numerics, 
physics, data) have to be identified, sized and solved (in 
turn) according to their (separate) contributions to a 
forecast. The goal, of course, is to prevent such problems 
from masking the effects of SEASAT data. Although the 
primary objective of this contract effort was to examine the 
effects of grid resolution on analysis and prediction using 
real data, its nature and scope became much broader as the 
work progressed. In addition to being able to document such 
grid effects, ODSI can also make a statement about the 
effects of: (1) variations in analysis procedures? and (2), 
variations in many discretionary procedures and parameters 
in prediction. Finally, this modeling experience has signifi- 
cantly increased the probability that valuable SEASAT data 
will not get "washed out" of modeled atmospheric structui'es? 
i.e., the computational viability of such data will be 
greatly improved as a consequence of tlie insight gained in 
this contract effort. 
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II. SCOPE OP W03SK 


The original Statement of Workf as contained in Article 
1 of JPL Contract No. 954668, provided that: 

• ODSI shall complete and test an analysis "scheme" 
(i.e., a set of computer programs tc perform 
objective analyses of certain atmospheric para- 
meters) needed to initialize four versions of an 
atmospheric prediction model. The design should 
permit the handling of both real and simulated 
observations . 


• ODSI shall develop and test tliree (out of four) 
versions of a primitive-equation prediction model. 

The grids to be used are: 63 x 63 x 5 layers? 

63 X 63 X 10 layers; 187 x 187 x 5 layers. The 
model should contain appropriate attributes 
(numerics and physics) for generating " opera tionally- 
competitive" forecasts of periods of one or two days. 


• ODSI shall exercise the analysis and prediction 
models with real data for two meteorological 
scenarios, the observations for which shall be 
provided by Fleet Numerical Weather Central. 


• ODSI shall verify, describe and document forecast 
model performance in both an objective (statistics, 
charts) and subjective (narrative) manner, empha- 
sizing those differences caused by variation in 
resolution. 


• ODSI shall perform such model development and 
testing on either the CYBER 175 (at PNWC) or the 
CDC 7600 (at NASA Ames) to the extent tnat computer 
resources permit. 
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This Statement of Work was subsequently modified to 
include : 

• The development and testing of a fourth version 
(187 X 187 X 10 layers) of a prediction model. 

. A task to conduct sensitivity tests with the 
coarse-mesh (63 x 63 x 5 layers) prediction 
model in order to determine and order the 
factors which might mask or impair the utility 
of SEASAT data. 

A. General Certification 

The terms and conditions of the Statement of Work (as 
amended) have been fulfilled completely. In fact, the scope 
of work actually performed goes well beyond the provisions 
of the statement, especially in those activities undertaken 
to develop/refine the outputs of the analysis and prediction 
models. The Scenario A data set was used several times, for 
example, to determine the cumulative/interactive impact of 
model changes on a sequence of twelve-hourly analyses and 
forecasts. Another example would be the extra efforts to 
introduce and to examine the effects of procedural changes 
in the models in order to improve the quality of model out- 
puts. Many of these improvements were incorporated in the 
final analysis production series. 


Vfith. respect to the forecast model runs, only eight 
prodxiction forecasts were needed to satisfy the axigmented 
work statement r since there are four model versions and two 
data sets. In practice ^ about 30-40 production forecasts 
were made. Most of these were coarse-mesh model forecasts, 
because not enough CDC 7600 time was available to perform 
extra runs with the 187 x 187 x 5 or 187 x 187 x 10 model 
versions. 

During the contract period, it became necessary to 
shift the programming/production efforts from one typo of 
computer to another (in order to accoimflodate GFE availa- 
bility) . Thus, some of the total effort was related to such 
system changes. 

B, Definitions and Cross-References 

Table II-l is the “x'oad map*' for review’ing Volxmie I. 

It contains descriptors, identifiers, and cross-reference 
information needed to facilitate the examination of study 
results. With respect to the figures, tables and charts 
which are referenced. Table II~1 is not intended to be 
exhaustive. 
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1- Analyses 

The analysis descriptors and selected chart refer- 
ences are located at the top of Table II-l. There are two 
•neteorological scenarios for which time i objective analyses 
were produced. For each scenario, twelve-hourly analyses 
were produced for each of two grids (63 x 63 and 187 x 187) . 
Production series analyses are interactive sequences of 
analyses (those which are coupled with first-guess fields 
produced by a 12- hour forecast from the previous analysis) . 
Others are "test series" ainalyses for a single date-time 
which are non- interactive analyses (the cumulative effects 
of the analysis-prediction sequence were not determined) . 

As indicated earlier, analyses were made using the 
187 X 187 grid, but these were not actually employed to 
initialize the fine-mesh forecasts for two important reasons: 
(1) the comparison of fine-mesh and coarse-mesh forecasts 
would be affected (in ways not easily assessed) by the 
presence of additional (small) scales in the fine-mesh 
analyses; and (2) there was not enough GDC 7600 computer 
time to run both types of fine-mesh forecasts (those using 
interpolated coarse-mesh initial values and those using fine- 
raesh initial values) . The fine-mesh forecast models were 
initialized, therefore, with interpolated coarse-mesh values 
to facilitate comparisons. 
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During tli<* analysis development effort, several 
test and production series analyses were improved upon and 
discarded* This included Production Series h and B, and 
Test Series A, One of the difficulties with early analysis 
versions was related to the procedure used to "assemble" each 
datum to the nearest grid point. Later versions allow each 
datum to influcsnce more tiian one grid point according to an 
area influence function. 

Selected ODSI analyses are shown in Section VIII. 
Charts VIII-1 through VIII-12 are key analyses in Scenario A, 
Charts VIII**13 through Vlll-24 are key analyses in Scenario 
B. These may be compared to FNWC analyses provided as Charts 
VIII-25 to VIII-SD, and VIII-31 to VIII-34, respectively. 

2. Forecasts 

Table II-l also provides model descriptors and run 
nunxbers for both the five-layer and ten-layer forecasts 
treated in this report. These are distributed by scenario. 

In addition, the table indicates the type of analysis model 
outputs used for forecast model initialisation, as well as 
the key refei'ences (tables, figures, and charts). 

The major runs for Scenario A are F18 (63 x 63 x 5) 
F22 (187 >: 187 x 5) and ^ (63 x 63 x 10) . Test Series B 
analysis outputs were used to initialise all three forecasts. 


(The 187 X 187 x 10 forecasts will be treated in Volume V 
of this Final Report.) Other runs in Scenario A are pro- 
vided for assessing the impact of (varying) analysis pro- 
cedures on short-range forecasts. Thus, P18 can be compared 
to F20 to assess the impact of analysis-forecast interaction. 

For Scenario B, Production Series C analyses were 
used throughout. The key runs are F19 (63 x 63 x 5) , F24 
(187 X 187 X 5), and T5 (63 x 63 x 10). 
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AND CROSS-REFERENCES 


SCENARIO A 

SCENARIO B j 

63 

187 X 187 

63 X 63 

187 X 187 

PRODUCTION 

TEST 

PRODUCTION 

PRODUCTION 

SERIES 

SERIES 

SERIES 

SERIES 

c 

B 

C 

c 


4 

13-15 

16 

- 

8 

17-19 

20 


“ 

99-102 

— 

PECHCV 

PEFHCV 

PECHCV 

PEFHCV 

F20 

F22 

Fl9 

F24 


4,43 

23,45 

24,46 

15 

12 

19 

18 

146-151 

131-133 

161-166 

158-160 

83-90 

63-66 

107-114 

103-106 

j 

PECUFV 

H/A 

PECHFV 

N/A 

T6 


T5 




22,44 


16 


20 


152-157 


167-172 


91-98 


115-122 










III. ANALYSIS MODEL DEVELOPMENT 


Objective analysis models were designed, developed and 
tested (on the PNWC CYBER 175) using real observational data. 
Each of the analysis programs may be executed on either the 
63 X 63 coarse-mesh grid (CMG) or the 187 x 187 fine-mesh 
grid (FMG) . These CMG and PMG polar-stereographic mappings 
of the Northern Hemisphere are true at 6 ON where the mesh 
distances are 381 and 127 kilometers, respectively, Elsev;here, 
the map factor 


m((|)) 


1 t sin 60 
1 + sin 4> 


is required to obtain true-earth distances. 

The parameters for analysis include: sea-surface 
temperature; twelve levels of wind, temperature and geopotential 
height (from 1000 to 100 millibars) ; and sea-level pressure. 
Stratospheric (50 to 10 millibars) temperatures and heights 
were generated using regression equations. Before objective 
analyses could be performed at the 950 MB and 900 MB levels, 
it was necessary to merge/check significant-level and mandatory- 
level radiosonde information for each analysis time. 

The primary objective of this modeling effort is to 
examine the effects of grid resolution on analysis and short- 
range weather prediction. A short discussion of and rationale 
for the model resolutions chosen will be presented in the 
following paragraphs. 
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The writers assume that the reader v;ill become familiar 
with the technical aspects of the analysis models as presented 
in Volume II of this Final Report. 

1. Vertical Resolution 

With respect to the lOOO-lQO MB mass structure analysis, 
FNWC partitions each vertical column into eight layers — 
each of which is assigned a value of static stability. In 
that model, the lowest layer (from 1000 - 775 MBS) is quite 
gross. The shallow stable layers (inversions) in the lower 
troposphere cannot be depicted. ODSI retained the upper seven 
layers but divided the 1000 - 775 MB layer into three sub- 
layers, each of which is assigned a value of static stability. 
The new layer interfaces occur at 925 MBS and 875 MBS. 

Thus, ODSI designed a ten-layer, twelve-level mass structure 
model. This led to the aforementioned requirement for 
inputs at 950 MBS and 900 MBS, 

2. Horizontal Resolution 

In the horizontal, ODSI selected the one-third mesh 
(187 K 187) grid for both analysis and prediction. On this 
grid, the true-earth mesh distance varies from 68 kilometers 
at the equator to 90 kilometers at 45N to 136 kilometers at 
the pole. It permits significantly better depiction of 
smaller scales in analysis. It dramatically reduces spatial 
truncation error in prediction. Indeed, for the shorter 
baroclinic waves the undermovement should be no greater than 
one or two percent. 
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A. 


Objective Analysis: Discussion 


The procedure (or scheme) used in objective analysis 
of atmospheric-state parameters is to invoke a set of rules 
to obtain a value at each discrete grid point by blending 
several types of information: (1) new observations (which 
are irregularly-spaced, non-representative of the scales 
being modeled, and contain as assortment of errors) ; (2) 
guess-field functional values and differential properties; 
and (3), control-field information (climatology). [If vari- 
ational techniques are employed, it is possible to impose 
dynamical constraints on the relationship between two or 
more parameters for analysis.) The relative impact of each 
type of information may be specified according to a set of 
weights. The objective of an atmospheric analysis is to 
describe/depict the "relevant** scales of motion, processes 
and effects consistent with problem characteristics, available 
resources (computer power, tiiwe, data) , and the needs of the 
users. In numerical weather prediction, one performs an 
analysis in order to pose the initial state. 

An objective analysis scheme tends to follow the same 
rules that a trained meteorologist follows as he produces a 
subjective (manual) analysis. An objective analysis, however, 
should be devoid of day-to-day, person-to-person variations 
in quality. Since the first-guess field embodies all of the 
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retained information of the parameter to be analyzed prior 
to receipt of current observations, both methods should and 
do exploit that resource. Since an analysis is primarily 
concerned with the locations and shapes of features (as 
refined through use of gradient and curvature information) , 
both methods may downgrade the importance of any single 
reported value. In some ways, the task is easier for a 
meteorologist. He can arbitrarily modify the procedure, 
smooth a given contour, or assign a measure of confidence to 
a report. These things caai be done in objective analyses, 
but not as easily. 

The grid resolution is extremely important in the 
objective analysis because it tends to define the range of 
spatial scales for treatment. Fine-resolution grids permit 
the depiction of more of the total atmospheric variation, 
but they also increase the accommodation of non-representative 
and/or erroneous reports. Pew, if any, surface observations 
are truly representative of the spatial scales of interest 
in NWP. The "local” content of observations tends to be 
aliased into the larger, permissible scales of motion. As 
the model resolution is changed, so must the grid-dependent 
procedviros and coefficients. The impact of each report on 
the analysis, for example, should depend more on the scales of 
features being represented than on the grid size. Thus, it 


becomes important to determine the influence functions for 


various observational classes and analysis types. Finally, 
an objective analysis should be able to introduce and distri- 
bute information into spatial scales which are computationally 
viable (in the prediction model) , and to enhance that viability 
through dynamic coupling and/or numerical consistency. 

Failing this, small-scale infonnation may get "washed out" 
in the adjustment period of a primitive-equation forecast 
model . 

B. A Scheme for Objective Analysis: Pattern 

Conservation Technique (PCT) 

ODSI has completed an objective analysis scheme, and 
has tested the scheme using a real data base assembled in 
conjunction with the Fleet Numerical Weather Central (PNWC) . 

All of the analysis programs use this Pattern-Conservation 
Technique (PCT) , which preserves specified differential 
properties of the first-guess field while fitting the obser- 
vations under a system of weights. ODSI selected this 
scheme because it exhibits many of the desirable character- 
istics of a good manual analysis procedure. 

The goal of a PCT scalar analysis is to blend the 
following information: the new data? the most recent past 
analysis or forecast value (the first guess) ; the gradients 
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of the first guess field in eight directions frora each grid 
point; and the Laplacian of the first guess field. The 
contribution of each such piece of information or term is 
specified by a set of weights. 

The desired fit of all analysis elements is realized 
fay niiniraising the sum of the squared deviations of the 
various characteristics of the analysis from their countex*- 
parts in the first guess. The minimization is accomplished 
with an elementary application of the calculus of vai'iations. 

Information is spread through space by the gradient and 
Laplacian terms. In a sea-level pressure analysis, for 
example, there are sometimes natural obstacles (mountain 
ridges, coastlines) beyond which a meteorologist would not 
allow a new observation to influence the analysis. This 
kind of constraint can be simulated in the objective analysis 
by reducing the influence of the gradients and Laplacian 
along the demarcation zone. 

Each PCT analysis "cycle" consists of three steps: (1) 
assemble the data at grid points using some specified influ- 
ence function; (2) solve an appropriate minimization equation; 
and (3) , re-evaluate the weight of each report. At least two 
cycles are needed to re-evaluate the weights. Indeed, an 
additional cycle (for a total of three) will often cause a 
rejected/dov'mgraded report to re-enter the analysis. 
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Vector applications of PCT may be regarded as extensions 
to the scalar procedures. In the objective wind analysis, for 
example, the vector application requires solution of two mini- 
mization equations (u-component and v-component equations) 
simultaneously. Also, the method preserves the vorticity and 
divergence of the first— guess wind field while fitting reported 
winds (under a system of weights) . 

The reader should review the contents of Sections I and II 
of Volume II which describe the scalar and vector applications 
of PCT in analysis. 

1. Solution Method for PCT Equations 

In the previous section, it was pointed out that 
the PCT scheme required solution of a "minimization equation" 
peculiar to the analysis type. Originally, the method of 
successive over-relaxation (SOR) was used. When tests were 
made to determine the impact of a single observation on the 
analysis, it was obvious that the SOR method was not distri- 
buting the information as would be expected from numerical 
and meteorological considerations. 

Several approaches were made to the solution of the 
minimization equation. The criteria by which the success of 
an approach was measured were of two kinds . One was the 
speed of convergence of iterative methods together with the 
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length of time until convergence. The second was the 
degx~ee to which the solution obtained agreed with meteoro- 
logical expectations in the three cases: using only one 
data report for a constant guess field? using only four data 
reports for a constant guess field (one report in each 
quadrant; and using a complete (typical) set of reports. 

The methods of solution which were explored were: 
Poisson fast solve method? successive over-relaxation by 
points? successive over-relaxation by lines? Jacobi iter- 
ations; Gauss iterations? block iterations? successive over- 



relaxation in 5-block scanning order; and successive over- 
relaxation scanning by spiral (the relaxation methods for 
the vectoi' minimisation equations required under-relaxation) . 
Further description of each method follows. 

The Poisson fast solve method uses the routine 
POIS 2 developed at NCAR. Although the minimisation equation 
is not a Poisson equation — or even elliptic — it is 
possible to extract a Poisson part of the equation and to 
attempt to solve it by iteration. More precisely, and in 


general, given the linear matrix equation Ax = B, one may 
write A = P - Q and the equation may be written Px = Qx + B. 
If the matrix P is amenable to some quick solution method, 
then one may attempt to find a solution via iteration: 





In tlae Poisson method, we choose P to he the Poisson matrix 
part of the equation. As is well known, a solution by iter- 
ation is possible for all initial guesses of x if and only 

-1 

if the eigenvalues of P Q are all of modulus less than 
unity. This was not the case for the minimisation equation, 
and iterative atten\pts at solutions diverged. 

The method of Jacobi iteration is a simultaneous 
relaxation method. In terms of the geneoral method above, 
we write our equation as Ax = 0 and A = (a. .) as 

^ / J 

(a. .6. j - a. .), where <S. . = [? This method failed 

1 ^ JL if J if J if] i i— J 


to converge. 

Gauss iteration is Jacobi iteration but using the 
updated values of the unkno^m as soon as they are available. 
This corresponds to splitting AasA= (L+I) +u where 
L is tlie lower triangular part of A, u is the upper trian- 
gular part, and I the main diagonal. Gauss iteration con- 
verged, but quite slowly. This had the undesirable effect 
of spreading data too far from the location of the report. 

The method of successive over-relaxation by points 
is an accelerated Gauss iteration method. In general, an 
interative method may possibly be accelerated by taking the 
splitting A = P - Q and rewriting it as A = (P + otQ) - 
(Q + ctQ) for some choice of real number ct H ~1- a 
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suitable choice of ct/ this method convei’ges quite ^ apidly# 
but nevertheless, allows data to spread in anomalous ways. 

Successive over-relaxation by lines is as above 
for points, but instead of updating the unknown as soon as 
each new x. is computed, update only after a line of x. ’s 

JL 

corresponding to a line of grid points is computed. This 
method failed to converge. 

The method of block iteration involves a splitting 
A = P - Q where P is block pentadiagonal. Since a block 
pentadiagonal matrix equation can be solved directly and 
rapidly, this method is suitable. It yielded the fastest 
convergence of all in the number of iterations, but took 
a longer execution time than the fastest method. In addition, 
there was a somewhat anomalous pattern in the spread of the 
data. 

The last two methods are both successive-over- 
relaxation but differ from the above in the order of the 
unknowns in the equation. This involves permuting rows of 
the matrix of the equation — but because of the splitting — 
this yields different effects. 

The 5-block scan involves grouping the unknowns 
into 5 groups and scanning them in order in each group, from 
group 1 to group 5, This method converged fairly rapidly — 
but not most rapidly — and yielded some anomalous data spread. 


K 
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The method of spiral scan was designed and tested. 
It requires scanning the grid from the center (pole) point 
in a counter-clockwise spiral outward to the edges. It 
converged fastest (of all) and yielded the least anomalous 
pattern of data spreading. 

Thus, the best result was obtained by SOR using 
a spiral scan with variable over-relaxation factor. Indeed, 
under-relaxation was found to converge more quickly in 'he 
simultaneous solution of the u- and v-component minimization 
equations in the objective wind analysis. This method has 
been implemented in all PCT codes in current use. 

2. Analysis Code Sequence and Computer 
Requirements 

The entire set of analysis codes outputs about 
fifty fields per horizontal grid per analysis time. The 
sequence of analysis codes commences with the 12-hour fore- 
cast outputs of the coarse-mesh PE Model (PECHCV) as init- 
ialized with analysis outputs for the time twelve hours 
previous to the current (analysis) time. Each program in 
the sequence is necessary either for a subsequent analysis 
code or for the forecast model, as listed below; 

a. Program to raerge/check radiosonde observations, 
and to interpolate the given profile to those pressure levels 
not already included in tlie observation; 
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b. Program to analyse the sea-surface tempera- 
ture (SST) ; 

c. Program to generate the first-guess sea- 
level pressure and upper-level temperature fields; 

d. Program to analyze the sea-level pressure; 

e. Program to analyze temperatures at twelve 
pressure levels from 1000 to 100 MBS; 

f. Program to generate the fir'st-guess upper- 
level geopotential-height fields (using outputs from d and 
e^ above) ; 

g. Program to analyze the geopotential-heights 

at twelve pressure levels from 1000 to 100 MBS, to "retrieve" 
hydrostatically-consistent virtual temperature pi'ofiles (if 
desired) , and to generate 50- , 30- and 10-MB temperatures 
and heights using regression equations; 

h. Program to analyze the reported winds 
(satellite cloud vectors, rawinsondes, aircraft reports, 
pilot balloons) at twelve pressure levels from 1000 to 
100 f4BS; 

i. Program which assembles the forecast model 
input fields. 

A substantial amount of CYBER 175 computer time 
was needed to develop these analysis codes. For the contract 
year, analysis model development and production runs consumed 
about 150 CP-hours (about 400 wall-clock hours) . 
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The program sequence shown above (exclusive of 
the 12-hour forecast) uses about six CP-minutes (twenty- 
three wall-clock minutes) on the coarse-mesh grid, and 
thirty-three CP-minutes (100 wall-clock minutes) on the 
fine-mesh grid per analysis ' time. There are nine analysis 
periods per scenario, and two scenarios. Each analysis 
requires a coarse-mesh, twelve-hour forecast — which 
takes seven CP-minutes and seventeen wall-clock minutes. 
This forecast is used for analyses on both grids. Thus, 
each analysis scenario takes about 400 CP-minutes or about 
1,400 wall-clock minutes to complete. There are two 
scenarios. For the first scenario (April 1976) , the entire 
scenario was run three times altogether, in order to imple- 
ment and test the cmtiulative effects of certain procedures. 


C. Data Influence Region 


The analysxs procedure xn use at the outset of this 
study permitted a report to be assembled in the analysxs 
at the nearest grid point. On the coarse-mesh (63 x 63) grid 
the procedure appears to be satisfactory in that it tended 
to accommodate the observations while preserving the smoothly 
varying differential properties of the first-guess field. 
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Once the fine-mesh (187 x 187) analysis codes were 
exercised, however, it became obvious that the procedure 
had to be altered in such a way as to cause an observation 
to impact on the analysis in the same geographical/raeteoro- 
logical scales — regardless of the grid size. As a report 
was introduced in the 187 x 187 analysis, the information 
was being distributed to the same number of gridpoints as 
in the 63 x 63 analysis, but the geographical distance/area 
was much less tlian in the 63 x 63 analysis. 

Procedures were implemented and tested in both the 
63 X 63 and 187 x 187 analysis codes which distributed 
observational information within a circular neighborhood 
in a manner which may be specified. (This action then 
pointed out the need to re-evaluate the weights of obser- 
vations in a similar manner.) For the 63 x 63 code, the 
implementation of this change was rather straightforward. 

The radius of the influence circle was chosen to be three, 
two, and one (coarse) grid lengths for the first, second and 
thix'd analysis cycles, respectively. For the 187 x 187 
code, the same geographical areas were influenced, but the 
number of grid points affected could be as large as 296, 

112, or 32, depending on the cycle number. 
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The Effects of Varying the Weights 


The constants which determine the relative importance 
of observations and the guess-field (and its differential 
properties) were examined. Tests showed the weights to be 
grid dependent . They also showed that the amount of com- 
puter time needed to arrive at a solution can vary markedly 
with changes in the weights. This is especially so in the 
case of the weight assigned to the Laplacian terms (in the 
pressure analysis) . Trial and error constitutes tlie only 
way to determine the optimum weights. There are no guar- 
antees, however, that the minimisation equation will even 
converge for some choices of weights. 

Tests were made to see if one can determine the combin- 
ation of weights which produces a valid, visually-pleasing 
analysis In a computationally-ef ficient way. The results of 
these tests were not particularly encouraging. Indeed, it 
appears to be the case that a substantial quality-control 
effort is required to prevent deterioration of the analyzed 
structures in data-poor areas. Two special cases can be 
discussed. First, it should be stated and understood that 
the primitive-equation forecast models do not predict the 
observed changes in the equatorial/subtropical latitudes 

ft 

with any degree of skill. Systematic, if only slight, 
errors are common in low latitudes where the initialization 
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problem is acute. Over an analysis-prsdiction-anaiysis 
scenario of several days in length, the gradual departures 
fron reality can become rather large in data-poor regions. 

If, after several days, an observation ^ received, it may be 
sufficiently different from the value (s) in the first-guess 
field as to be doxmgraded or discarded. The opposite situ- 
ation is also troublesome. That is to say, an analysis scheme 
may accept a report xdiich, because of its difference from the 
first-guess values, may lead to a perturbation in tlxe analysis 
(it could be real, but it is probably spurious) . The primi- 
tive-equation forecast model x^ill, more than likely, amplify 
such a perturbation into quite an intense, unreal circulation. 
The problem, of course, is to develop techniques for distri- 
buting observed information into appropriate ranges of scale 
(wholesale updating of guess-field values) rather than to 
allox\r a report to create a spurious oirculation/system. If 
such techniques cannot be developed, it xvould be better not 
to introduce isolated data at all in sparse-data regions 
which are characterised by X\feaJ; horizontal gradients (generally) 
and large horizontal scales. 



E. 


Some Practical Solutions to Analysis Problems 


As in any undertaking of this scope and complexity ^ 
many problems were encountered, and either solved or cir- 
cumvented* At first, the problems tended to be related to 
the implementation of model procedures in terms of correct/ 
efficient code. But, once these difficulties were overcome, 
the problems centered on the effects of sparse and/or bad 
data (especially in the tropics) and on the procedures for 
distributing information in isolated reports (spatially) 
into larger meteorological scales. On the other hand, no 
amount of effort can solve the fundamental problem in 
objective analysis: how do you handle an unsupported obser- 
vation of an unexpected change/ event? Another important 
problem arises from any procedure which relies too heavily 
upon either the reports or upon the guess field. Even 
though it can be argued that the guess field represents 
the cumulative knowledge of the atmosphere prior to receipt 
of the latest observations, the guess field often departs 
from reality in data-poor regions so much so that perfectly 
good observations will be downgraded or rejected altogether. 

The use of (satellite) cloud photographs, therefore, is vital 
to the required quality control effort in operational analysis . 
Even so, the design (and redesign) of data toss-out criteria 
(with provision for manual override) represents an important 
and open-ended aspect of analysis modeling. 
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Problems in Low Latitudes 


During the early stages of the analysis development 
effort, analysis outputs contained unrealistic, small-scale 
features in low latitudes, generally, and in the lateral 
boundary areas, specifically. SIRS soundings caused the 
general problem, because the actual horizontal gradients in 
the tropical regions tend to be small compared to ficticious 
gradients caused by errors in the SIRS measurements. Once the 
small-scale features get into the analysed structure, the 
primitive-equation model has a tendency to produce (forecast) 
structures which are even less realistic. Since tlie analysis 
sequence (over several days) is maintained with 12-hour fore- 
casts, each analysis can become increasingly unrealistic. 

Charts VIII-49 and VlII-50 illustrate the loxir- 
latitude analysis problem. Chart VIIl-49 is a 700 MB tempera- 
ture analysis for 12002, 15 April 197G. It contains numerous 
satellite temper'ature soundings. Chart VIIl-50 is the cortes- 
ponding 700 MB height analysis containing SIRS height obser- 
vations, An early version of the analysis procedure was used 
for these outputs. The non-meteoro logical character of these 
outputs is clearly indicated in the ocean regions south of 
(about) 30° North, 


111-18 


Ghaarfes and VllI-52 wave also generated 

vising the sanio (early) versions of tho analysis models. 

Chart vm-51 contains tho sea-lovel pressure analysis for 
ISQO^i 22 hpril 1976. It is quite noisy in appearance, 
ospoeially in low latitudes. Chart VllI-52 contains the 
sea-lev 'l pressure analysis twenty- four hours later (1200S, 

22 hpril) . Reeause of the intoraetion with the forecast 
model (described above) , this pressure pattern has become 
quite unrealistic generally, but totally unacceptable on 
the lewer/upper boundaries. This illustrates how seemingly 
innocuous errors (shown in Chart VIII-51) can amplify after 
only (subsequent) 12-hour forecasts and analyses (each) . 

Several remedial procedures wore implemented and 
tested in the analyses. These included: (1) Fourier filtering; 
(2) Shapiro filterinvi; (3) Laplacian-typo smoothing with 
latitudinal variation; and (4), removal of SIRS reports 
south of a proscribed latitude (about 20N) . The first two 
measures were discarded. Some remedial steps wore also 
taken in the forecast model. These included; (1) use of a 
tendency truncator which does not permit forecast changes 
south of 5N, but which leaves the tendencies unaltered north 
of ISN (with linear variation between these limits); (2) 
application of a low-pass filter to model outputs (with no 
feedback to the integrations) ; (3) uso of thermal/momentura 
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diffusion in the integrations; and (4) , use of latitudin- 
ally-dependent Laplacian-type smoothing (between 5N and 30N) 
of model outputs. 

Chart VIII-2 illustrates the cumulative benefits 
of these remedial actions. This 1200Z, 23 April PS analysis 
should be compared to Chart VIII-52. On the basis of many 
model tests ^ it was concluded that the entire 63 x 63 analysis- 
prediction scenario could be carried out without further modi- 
fications to the analysis procedui-e. 

2 . Analyzed vs. Retrieved Temperatures 

The original analysis technique called for modeling 
the 1000-100 MB mass structure in terms of eight (or ten) 
static-stability layers. Once the heights, Z(p), were 
determined for the column, the solution of a matrix trans- 
formation equation yielded a corresponding set of vii'tual 
* 

temperatures, T (p) , which could be integrated hydrostatically 
to produce the heights, Z{p). These so-called "retrieved" 
temperatures tend to diffei", however, from analyzed tempera- 
tures at any given pressure level. 

Chart VIII -5 3 shows the retrieved temperatures at 
950 MBS for 1200Z, 22 April 1976. In Chart VIII-54, we show 
the 950 MB analyzed temperatures. These shox^; that the 
retrieved temperatures are not only noisy in appearance, but 
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they also lose much of the structure normally found in care- 
fully-drawn subjective analyses. Analysed temperatures were 
used thereafter to initialize the prediction model with no 
apparent negative effects. 
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Data Base Preparation 


The analysis program sequence described in Section 
III.B.2 was executed using real observations in order to 
prepare the data base needed to Initialize the several 
versions of the prediction model. The CYBER 175 at PNtfC was 
utilized for all aspects of observation processing and 
objective analysis. 

There are tvro meteorological scenarios of four days 
duration each. Scenario A covers the period from 1200Z, 21 
April 1976 through 12002, 25 April 1976. Scenario B covers 
the period from 12002, 19 May 1976 through 12002, 23 May 
1976. The analyses are performed at twelve-hourly intervals. 
Thus, th^re are nine analyses in each four-day scenario. As 
noted earlier. Scenario A was the subject of special investi- 
gations leading to improved procedures in analysis. The 
analyses were performed on both the 63 x 63 coarse-raesh and 
the 187 X 187 fine-mesh grids. 

Table II-l indicated that the Test Series B analyses 
were used to initialize the forecasts in Scenario A. These 
are non-interactive analyses. Specifically, analyses were 
made for individual date-times using an area-influence pro- 
cedure without allowing the forecast model to propagate such 
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effects to a subsequent analysis. In contrast, the Production 
Series C analyses were fully interactive. (By examining the 
P20 and PIS forecasts, one can assess the impact of such 
interactive procedures on a short-range forecast.) 

The scenarios were selected for periods during which 
interesting (and large) changes were taking place. All 
available conventional and satellite data at PNWC, together 
witli corresponding PNT^C analyses and forecasts, were 
collected with the help of PNWC personnel. The data counts 
for 12003, 22 April 1976 are representative. Included are; 

(1) 4,706 sea-level pressure reports; (2) 442 sea-sux'face 
temperatures; and (3) , upper-air tempex-atui'e , height and 
wind repox“ts in numbex's shown in the table below; 


Number of Reports 


Pressure 

Tempex'ature 

Height 

Wind 

1000 

402 

476 

266 

950 

530 

496 

326 

900 

562 

504 

691 

850 

637 

647 

MSG 

700 

651 

654 

551 

500 

644 

649 

550 

400 

637 

640 

557 

300 

631 

632 

604 

250 

MSG 

MSG 

899 

200 

593 

595 

976 

150 

566 

567 

476 

100 

489 

491 

470 


\ 
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At this date- time, there were about 500 radiosonde/rawinsonde 
soundings, but values at specific levels were sometimes missing 
or garbled. CouBfts in excess of 500 generally indicate the 
presence of SIRS soundings/winds. Satellite wind vectors are 
available at 900 J4BS, and in large numbers at/near the 200-MB 
level. 


G. Charts and Tables 

Recall that the analysis models output approximately 
fifty distributions/fields at each analysis time on each 
horizontal grid. Thus, the rrharts provided in this Final 
Report represent only a small fi action of the total. In 
addition to presenting selected analyses, analysis actual 
change charts (the difference between analyses at tv/o times) 
and analysis difference charts (the difference betxi?een analyses 
on two grids) will be presented, as appropriate. Statistical 
measures of such changes/differences will also be provided. 

1. Scenario A Discussion 

Although this scenario commences at 1200Z, 21 April 
1976, the forecasts were initialized with analysis model out- 
puts for 1200Z, 22 April 1976. 

Charts VIII-1 through VIII-3 contain the sea-level 
pressure (PS) analyses for 22 April through 24 April, respect- 
ively. These were produced on the 63 x 63 grid. During this 


period, the atmosphere Tinderwent some rather large changes. 
Chart VIII-59, for example, shows the PS actual changes 
{using a 4-MB contour interval) during the first day. There 
are several change centers with magnitudes greater than 
taventy millifcars (about 1 MB per hour) . Charts VIII-5 through 
VIlI-7 contain the 500 MB height analyses for the 22-24 April 
period. The main feature is the blocking high southeast of 
Iceland, which leads to height falls over southern Europe 
(do\«istream of the block) . In contrast, the Pacific region 
undergoes some rather large adjustments in the open-wave 
pattern. The Pacific Ocean atmospheric changes are clearly 
more interesting than those in the Atlantic . For the most 
pctrt, they reflect system intensity changes. In contrast, 
the intense low pressure system in the Atlantic (near 45N/ 

40W) changes intensity and location only slightly. In fact, 
the Atlantic and European sectors are comparatively inactive 
during this period. Chart VIII-60 contains the corresponding 
changes in the 500 MB heights. This shows a maximum change 
of -360 meters over China, with several change centers in 
excess of + 180 meters. The 48-hour changes are shown in 
Charts VIII-61 and VIII-62. 

Table VII-11 contains some statistical measures of 
the atmospheric changes taking place in this 48-hour period. 
The RMS pressure difference increases from 4.81 MBS to 
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6.01 MBS after one and t\>ro days^ respectively. The RMS 500 MB 
height difference increases from 54.7 meters to 67.7 meters 
after one and two days, respectively. 

In the evaluation of prediction model outputs, 
emphasis will be given to the system over China and to the 
system in the Gulf of Alaska. The former deepened 19 MBS in 
the first day, and the latter deepened 23 MBS in the first 
day. The Atlantic low may be used as a standard since 
it did not change in intensity during the period. 

Chart VI I I -9 contains the 500 MB temperature 
analysis for 1200Z, 22 April 1976. Note the SIRS reports 
over the oceans north of 2 ON. 

Chart VIII-IO contains the 500 MB wind analysis 
for 1200Z, 22 April 1976. Reports are shown as heavy barbs. 
Chart VIII-12 contains the 900 MB wind analysis for the same 
time. This contains many satellite wind vectors. 

The reader may compare the aforementioned PS and 
Z500 analysis outputs to FNWC analyses. See Charts VIII-25 
and VIII-26 for the PS charts for 22 and 23 April, respect- 
ively. See Charts V3 11-27 and VIII-28 for the corresponding 
500 MB height analyses. 


2 . 


Scenario B Discussion 


Although this scenario commences at 1200Z, 19 May 
1976, the forecasts were initialized with analysis model 
outputs for 1200Z, 20 May 1976. 

Charts VI 11-13 through VI 11-15 contain the sea- 
level pressure (PS) analyses for 20 May through 22 May 1976, 
respectively. These were produced on the 63 x 63 grid. 

Charts VIII-99 and VIII-101 contain the PS actual changes for 
the first day and total period, respectively. The reader 
will note that a typhoon (present near 12N/146E) accounts for 
the rise-fall pair of centers on Chart VIII-99. 

Charts VIII-17 through VIII-19 contain the 500 MB 
height analyses for 20-22 May. Charts VIII-100 and VIII-102 
contain the Z500 actual changes (60-meter contours) for the 
first day and total period, respectively. For the most part, 
the actual changes for Scenario B are smaller tlian those 
indicated in Scenario A, but they are still large enough to 
be of use in this study. In Chart VIII-102, for example, two 
Z500 change centers exceed 300 meters in magnitude. 

As Charts VIII-13 and VIII-14 are examined, one 
finds that the intensity of centers did not change dramati- 
cally. One secondary center between Iceland and England 
deepened from 1004 to 996 MBS. The deep low over China 




III-27 




remained constant at 988 MBS. The low center east of 
Kamchatka and the low center in the Gulf of Alaska each 
filled about four millibars. In the central Atlantic (near 
40N/40W) , the low deepened a little from 1009 to 1006 MBS. 

Thus,, the displacement of systems was a significant consider- 
ation in this scenario. 

Selected additional parameters are also provided 
at the (base) time of 1200, 20 May 1976. Chart VIII-21 
contains a 500 MB temperature analysis. Chart VIII-22 
contains the corresponding wind analysis. Both were per- 
formed on the 63 x 63 grid. 

The charts for Scenario B may oe compared to FNWC 
analyses. Charts VIII-31 and VIII-32 contain the PS analyses 
for 20 and 21 May, respectively. Charts VIII-33 and VIII-34 
contain the Z500 analyses for these same dates. As would be 
expected, the correspondence between ODSI and FNWC analyses 
is quite good, especially in regions with adequate obser- 
vational counts. 

Table VII-17 contains some statistical measures of 
the actual changes at sea level and 500 MBS in Scenario B. 

The RMS pressure change amounts to 3.34 MBS in the first day 
(20-21 May), and increases slightly to 4.35 MBS for the two-day 
period (20-22 May) . These are somewhat smaller than the RMS 
changes (of 4.81 and 6.01 I4BS) in Scenario A, 
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H, Effects of Grid Resolution 

Selected analysis model outputs using the 187 s 187 grid 
are provided in order to facilitate comparison with the 63 x 
63 outputs. 

1. Scenario A 

Chart Vlll-4 contains the 187 x 18? sea-lovel 
pressure analysis for 1200S, 22 April 1976. It may bo oom- 
parod to Chart In fact, the (small) differences 

between these two analyses are contained in Chart VTII-65. 

(The contour interval is 1 MB.) 

Chart VIII”8 contains the 187 x 187 SOO MB l^eiaht 
analysis for 12003, 22 April 1976. Compare this with the 
63 X 63 version shown in Chart VII 1-5, The ditforencos 
betv\'oen these analyses are shown on Chart VIIT-5^. (The 
contour interval is 10 meters.) Here, one starts to see 
differences of 10-30 meters in places. Note that in deuso- 
data areas, the differences between the 187 x 187 and 63 x 63 
outputs occur ix\ small scales (as expected, bocavise of the 
ability of the PMG to accommodate those reports better) . In 
sparse-data regions, the differences appear to be in larger 
(metoorologicai) scales. This is due to the fact that; (1) 
filtering action on the 63 x 63 affects a bro.adev geographical 
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region than on a 187 x 187 grid? and (2) , the percent of the 
difference (between the reported value and the first-guess 
value) which gets applied to a grid point will tend to be 
larger for the 187 x 187 grid because of proximity (alone) . 

Chart VIII-11 contains the 500 MB temperature 
analysis for 12002, 22 April 1976 for the 187 x 187 grid. 
Compare this to Chart VIII-9. 

2. Scenario B 

Chart VIlI-16 contains the 187 x 187 sea-level 
pressure analysis for 1200Z, 20 May 1976. Compare this output 
to Chart VlII-13 containing the 63 x 63 output. The differ- 
ences between these two analyses are shown on Chart VIII-57. 
There are two small regions of 2 MB differences (central 
Pacific and Ethiopia). Otherwise, tlie differences are minimal. 

Chart VIII-20 contains the 187 x 187 500 MB height 
analysis for 12002, 20 May. Compare this to Chart VIII-17. 

The differences between these outputs are shown in Chart 
VIII-58. Once again, the differences are in very small scales 
in dense-data regions. In sparse-data areas (such as the 
central Pacific), the differences lie in larger spatial scales. 
The 90-meter difference near 45N/170W is the result of a 
single report (which produced a minor trough) in the 187 x 187 
analysis. 
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Chart VIII-23 contains the 187 x 187 500 MB tempera- 
ture analysis for 20 May. Compare this to Chart VIII-21. 

Chart VIII-24 contains the 187 x 187 500 i4B wind 
vector analysis for 20 May. . Compare this to Chart VIII-22. 

3. Sea-Surface Temperature (S5T) 

Eight sea-surface temperature analyses are pro- 
vided. Charts VIII-123 through VIII-126 pertain to Scenario 

A. Charts VIII-123 is the 187 x 187 SST analysis, which 
should be compared to Chart VIII-124 containing the 63 x 63 
SST analysis for the same date-time. The ability of the PMG 
to resolve strong gradients is extremely good. 

Charts VIII-127 through VIII-130 pertain to Scenario 

B. In this sequence, compare the 187 x 187 analysis in Chart 
VIII-127 with the 63 x 63 analysis in Chart VIII-128. 
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IV. PREDICTION MODEL DEVELOPifflNT 

The primary objective of this modeling effort was to 
examine the impact of grid resolution on short-range weather 
prediction (with real data) . In so doing, ODSI has prepared 
a more effective/efficient model context for assessing the 
impact of SEASAT data on weather prediction. The models 
possess many of the characteristics of existing operational 
models at Fleet Numerical Weather Central, Monterey. This 
will facilitate the applicability and transferability of study 
results to that computational environment. This is particularly 
true with respect to the procedures in analysis codes and to the 
physics in prediction codes. With respect to data management 
and computational procedures in the prediction models, ODSI has 
implemented alternative approaches in order to be able to exe- 
cute high-resolution versions on available computer systems. 

These tests have pointed out the importance of "engin- 
eering” in prediction model development. Computational devices 
(smoothers and filters, tendency truncators, mountain shaping, 
diffusion operators) are needed to control noise in the solu- 
tions. But, the repetitive use of such devices can/does 
produce effects which mask (or overwhelm) the effects of 
physical processes being modeled. Thus, the impact of each 
discretionary procedure needs to be ordered, understood and 
exploited. 
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Observations are deficient in both number and distri- 
bution. Ocean-area data coverage is a critical issue in 
objective analysis. Satellite measurements are intended to 
satisfy ocean-area data requirements eventually, but model 
development can/must proceed in the interim. When the full 
potential of this data resource is realized, there should be 
significant improvement both in the specification of the 
initial state for prediction and in our understanding of the 
underlying physics. The point, of course, is that special 
care must be taken to ensure the computational viability of 
observed information. 

A. Model Description 

A complete description of the various forecast model 
versions (and programs) is provided in Volume III of this 
Final Report. For reader convenience, however, some important 
attributes and procedures will be described briefly in the 
following paragraphs. 

The model requires a computer solution, using numerical 
methods, to a set of non-linear, partial difference equations 
at each grid point in the three-dimensional domain each inte- 
gration time step. The time step is 4 minutes (12 minutes) 
in the 187 x 187 (63 x 63) model versions. The horizontal 
domain is a polar-stereographic mapping of the Northern 
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Hemisphere. Terrain-following "sigma" surfaces are used to 
map variables vertically. PECHFV has ten sigma surfaces 
(a = 0.95, 0.85, ..., 0.05). PEFHCV and PECHCV have five 
sigma surfaces (a => 0.9, 0.7, ..., 0.1). 

The set of equations contains mathematical terms which 
represent relevant physical/dynamical processes and effects. 
The equations are written in flux form to facilitate the use 
of a second-order differencing scheme with special conser- 
vation properties. For temporal differencing, each twelve 
hours of integrations is initiated with an Euler-backward 
step and carried forward with centered steps (leapfrog) . 
Consecutive solutions are time-filtered using a procedure 
devised by Robert in order to minimise solution separation 
while selectively dissipating high-frequency computational 
modes. The pressure-gradient force terms in the momentum 
equations are time-averaged to permit a slightly larger inte- 
gration time step than would be acceptable ordinarily. 

Some of the more important physical/dynamical processes 
represented in the model ai*e: evaporation and condensation 
(cyclone-scale and cumulus-scale) sensible boat exchanges 
(air— ocean; air-ground) ; and solar and terrestrial radiation. 
The cumulus parameterization is based on a 1968 Araka\i?a 
algorithni, and redistributes heat and moisture in each column 
if certain (conditional-instability) criteria are met. Dry- 
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convective adjustments are invoked as necessary to preclude 
hydrostatic instability. Both mountain and frictional effects 
lend realism to flows. The quadratic friction term vanishes 
at sigma =0.8 (0.9) level in the 5-layer (10-layer) model. 
Mountain gradients cannot exceed 2,000 meters per meshlength. 

A land-sea-ice discriminator is employed, as appropriate. 

Terms representing the three-dimensional transport and 
vergence of heat, moisture and momentum are present, of 
course. Adiabatic temperature changes are modeled. 

The diabatic heating rates are calculated (each hour) 
in a manner similar to that used at UCLA, as outlined by 
Langlois and Kwok (1969) . Moisture is carried out at the 
lowest three (six) levels in the 5-level (10-level) models. 
Initial moistures are parameterized in terms of the patterns 
of geostrophic relative vorticity (the flow patterns) . 

Heating rates are calculated for two gross layers (between 
sigma = 1.0 and 0.6, and between 0,6 and 0.2). A gross 
cloud is modeled after Smagorinsky (1960) , and used to 
apportion the clear-sky and cloudy-sky rates. The solar 
calculations follow Joseph (1966) . Infra-red fluxes are 
based on Mintz and Arakawa, after approximate fluxes are 
calculated using the approach of Danard (1969) . 

The model is initialized using the objective analysis 
model outputs described earlier. Recall that vertical 
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interpolation is needed to initialize sigma-surface points 
using pressure-surface information. Winds come from an 
objective analysis in which the first-guess values are non- 
divergent geostrophic winds. 

Both explicit and implicit smoothing/filtering are 
present in this model. Without thermal and momentum diffusion 
terms, it would be difficult to control noise {hence, 
stability) . A non-linear pressure smoother is also employed. 

In addition, a tendency truncator is used to eliminate fore- 
cast changes south of 5° latitude, and to gradually increase 
(with latitude) the permissible amount of forecast change up 
to 20° latitude (above which all forecast changes are allowed) . 
[Such limiters are motivated and justified by problems with 
the data, (artificial) boundary conditions, and initialization 
procedures in tropical latitudes.] 

This class of model has been in operational use by FNWC 
since 1S70. It produces forecasts with skill (compared to a 
persistence forecast) for periods up to several days. Short 
baroclinic waves tend to be unde3nnoved by 10-20 percent on the 
63 X 63 grid, but only 1-4 percent on the 187 x 187 grid. 
Precipitation tends to be modeled well, especially in maritime 
cyclones. Indeed, the model handles cyclogenesis quite well, 
especially in the western Pacific Ocean. 



B. Special Characteristics 

Although a complete model description is provided in 
Volume III of the Final Report r brief explanations of special 
model characteristics will bo pr 'vided in this section. 

1. Geos tr ophi c 50 -MB Winds 

In both PECHCV and PEFHCV, the sigma = 0.1 surface 
is often above 100 MBS, especially over high terrain where it 
can reach to 50 JIBS. With PECHFV, the uppermost computational 
surface lies in the 50- to 25-MB stratum. Since the wind 
analysis model provides winds up to 100 MBS only^ one can 
either extrapolate from available levels or use a non-divergent 
geostrophic wind above 100 MBS. Tests showed the latter pro- 
cedure to be quite satisfactory {and an improvement over 
extrapolation) . 

2. Moisture Initialization Algorithm 

As pointed out earlier, the initial relative humidity 
(RH) values are parameterized using geostrophic relative vorti- 
city values from the flow patterns. This has proved to be 
superior to any so-called moisture analysis (which has few if 
any reports over the oceans) , or to climatology values. 
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Special Processing of Outputs 


Early tests showed that noise and/or unrealistic 
structures tended to build up in tropical regions during the 
course of the analysis-prediction cycle. 

The remedial procedures included restoring the 
initial boundary conditions for output distributions, with 
Laplacian smoothing (latitudinally dependent) south of 30° 
latitude. (In an operational context, this procedure would 
have to be replaced with one involving a return to climatology 
in analysis in data-sparse regions.) 


4. Special Processing of Inputs to Fine-Mesh 
(187 X 187) Models 

Because of the large data base with the 187 x 187 
grid (about one million words) , it is more efficient to sort 
the input (analysis) data into vertical sections (I,K) for 
which pressure is the vertical coordinate (same as in the 
analysis codes) prior to performing the vertical interpolation 
to sigma surfaces. A separate program named PEPREP performs 
this task and generates an appropriate file for use in the 
initialization overlay of the PE model. PEPREP uses both 
rotating mass storage and either ECS or LCK (about 250,000 
words) for temporary storage. 


The 187 X 187 model may be initialized with 
187 X 187 analysis outputs, or with 63 x 63 analysis outputs 
which have been interpolated to the 187 x 187 grid. 

5. Lateral Boundary Conditions and Tendency 
Truncation 

Forecast tendencies for each integration time step 
are reduced/eliminated south of 20N according to an array of 
coefficients. South of 5N the tendencies are set to zero 
(producing a persistence forecast) - Between 5N and 2 ON, the 
coefficients vary according to the latitude 0 as follows: 

P (sin 0 - sin 5°) ^ 

|^(sin 20'=’ ~ sin 5®) 

Within each integration step, however, the boundary walls 
are assumed to be rigid, insulated and slippery. No boundary 
fluxes are allowed. Tangential flow is permitted. 


6. Pressure-Gradient Force (PGF) Time-Averaging 
in Momentum Equations 

To allow for a slightly lengthened integration time 
step, the pressure-gradient force (PGF) terms in the momentum 
equations are time-averaged. The procedure requires that the 
integrations be re-ordered as follows: (1) integrate all of 
the equations except the momentum equations [in order to 



produce new (T+1) geopotentials]? (2) time-average the PGF 
terms, using 

= (1-a) PGP^ + I (PGr*^”^ + PGp'^'’^^)? 

and (3) , integrate the momentum equations containing the 
time-averaged PGP terms, 

7. Temporal Filtering of Thermodynamic and 
Moisture Equations 

In the temperature and moisture equations, consecu- 
tive solutions are time-filtered using a method devised by 
Robert (1966) . The method is used to minimise solution 
separation, as well as to damp undesirable short-wave com- 
ponents. It is not applied to the momentum equations. 

® • Non-Linear Pressure (Tendency Equation) 
Smoothing 

Oliger and Wellck (1970) describe a non-linear 
pressure- smoother which we apply in the pressure- tendency 
equation each time step. Since the operator is actually 
applied to stencils of sea-level pressure values about each 
point, the terrain pressure must first be reduced to sea 
level, then recovered from sea-level values, using a specified 
algoritlmi. 


9. 


Horizontal Diffusion 


Temperature diffusion, and moment- im diffusion, 

Dj^, were added to the model equations to redistribute/attenuate 
high-frequency components. The former (D^) varies in magnitude 
according to the diffusion coefficient, K, given by 

f 10*^ on the CMG 

^ = \ 5 

L 4x10 on the FMG 

The latter (D^^) varies with both latitude and model level, in 
addition to the aforementioned grid-dependent coefficient. 

10. Temporal Differencing 

Centered (leapfrog) differencing is used, except 
for the step wliich initiates a new integration period 
(normally every twelve forecast hours) ^ For this step, we 
use an Euler-backward step (which is selectively dissipative) . 

In the 187 x 187 model versions, a provision has been made to 
employ an Euler-backward step at other specified times, as well. 

Due to the interaction of the various computational 
devices being used, the integration is fifteen minutes in the 
63 X 63 models, and four minutes in the 187 x 187 models. We 
refer to: (1) PGP time- averaging; (2) Robert time-filtering; 

(3) non— linear pressure smciothing; and C4) lateral diffusion. 
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All of the analysis and coarse-mesh (63 x 63) prediction 
modeling was performed on the CYBER 175 at Fleet Numerical 
Weather Central ^ Monterey. In this contract period, this 
amounted to about 200 CP~hours and about 600 wall-clock hours. 
Table IV-1 shows the CYBER 175 utilization for typical 24-hour 
forecasts using either the 63 x 63 x 5- level or 63 x 63 x 10- 
level prediction models. 

The fine-mesh (187 x 187; prediction modeling was per- 
formed on the CDC 7600 at N.ASA Ames. Table IV-2 shows the time 
it takes to execute these high-resolution model versions. 

These tables indicate that it takes about 12-13 times longer 
to execute the one- third mesh version on the 7600 than it doe? 
to execute a standard-mesh model on the 175, Since the 7600 
is about twice as fast as the 175 for such codes, this 
suggests that the fine-mesh problem is a 25X problem. Theo- 
retically, it is a 27X problem. Thus, some economies of scale 
were realised. 

The lack of 7600 time prevented ODSI from executing the 
fine-mesh models on more data sets or for longer forecast 
periods. 

During the contract period, it became necessary to shift 
the programmi g/production efforts from one type of computer 
to another (in order to accommodate GFE at'ailability) . Thus, 
some of the total effort was related to such system changes. 
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TABLE IV-1: CYBER 175 UTILIZATION FOR 

COARSE-MESH PREDICTION. 



NOTES: 

1. Used 12-minute times teps. 

2. FKWC provided 200 CP-hour s (about 600 wall-clock hours) 
in this contract year. 




STABLE IV-2: 


CDC 7600 UTILISATION FOR 
FINE-i>tESH PREDICTION. 


MODEL/ ITEM 

CDC 7600 
EXECUTION TIMES 

WALL 

CP 

187 X 187 K 10 PEM (PEPHFV) 



— Initialisation 

44 Sec. 

13 Sec. 

— Each 6- hour Integration 

3.5 Hr. 

1 . 2 Hr . 

— Each 6-hour Output 

14 Min. 

1.5 Min. 

— Total 24-hour Forecast 

15 Hr. 

4.8 Hr. 

187 X 187 X 5 PEM (PEFHCV) 



— Initialisation 

\ 

23 .Sec. 

I 8 Sec, i 

1 1 

— Each 6-hour Intecrration 

! 67 Min. 

39 Min. ‘ 

1 

— Each 6-hour Outptit 

■ 9 Min . 

1.3 Min. 

— Total 2‘^-hour Forecast 

5.2 Hr. 

2.7 Hr, 


NOTES: 

1. Used 4-minuto timostops. 

2. Runs wore made at NASA Ames, 


V. PREDICTION MODEL TEST RESULTS 


Table II-l provides the model descriptors, run numbers 
and selected chart, table and figure references for those 
prediction model results to be presented in this Pinal Report. 
The coarse-mesh (63 x 63 x -5 or 63 x 63 x 10) forecasts are 
all two-day forecasts. The fine-mesh (187 x 187 x 5) fore- 
casts are one-day forecasts. There are three major runs in 
each of two scenarios. In Scenario A, all forecasts were 
initialized with Test Series E analyses. The major runs for 
Scenario A are: F18 (63 x 63 x 5) , F22 (187 x 187 x 5) and 
T4 (63 X 63 X 10) . Runs F20 (63 x 63 x 5) and (63 x 63 x 10) 

were initialized with Production Series C analyses, and may 
be compared to the major runs, as appropriate, to determine the 
effects of interactive analyses. In Scenario B, all forecasts 
were initialized with Production Series C analyses. The major 
runs for Scenario B are: F19 (63 x 63 x 5) , F24 (187 x 187 x 5) 
and ^ (63 x 63 x 10) . 

Several classes of exhibits will be provided to aid in 
the examination and evaluation of these results: (1) forecast 
charts; (2) forecast change charts (forecast minus starting 
analysis) ; (3) forecast error charts (forecast minus verifi- 
cation analysis) ; (4) forecast difference charts (ten-layer 
forecast minus five- layer forecast and fine-mesh forecast minus 
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coarse-mesh. forecast) ; (5) starting/verification analysis 
charts; (6) actual change charts (to be compared, to either 
forecast error charts or forecast change charts) ; (7) sta- 
tistical information describing model performance and model 
behavior (energetics) ; and (8) graphical information describing 
model behavior. 

A. Scenario A; Forecasts, Forecast Changes, Forecast 
Errors 

Several types of comparisons may be made in this scen- 
ario. The forecast charts may be verified against the 
appropriate verification analyses. Forecasts from a five- 
layer model (for a given initial analysis type) may be 
compared to forecasts from a ten-layer model. Forecasts 
(from a given model version) using one type of initial 
analysis may be compared to forecasts using a different type 
of initial analysis. Ideally , the forecast changes and 
actual changes should correspond quite well. Also, the 
pattern of forecast errors should be less dense/smaller in 
size than the pattern of actual changes. Such qualitative 
assessments utay be augmented by detailed comparisons of 
statistical measures presented herein. 
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1 . 


Run F22 (PEFHCV) 


Charts VIII— 131 through VIII~133 contain the 
24-hour sea“level pressxire (PS) , 500-MB height (Z500) and 
500-MB temperature (T500) forecasts, respectively, using the 
187 X 187 X 5 model version and lest Series B analyses. 

Charts VIII-2 and VIII-6 may be used to verify the PS and 
Z500 forecasts, respectively. The forecasts may be compared 
also to Charts VIII-134 through VIII-136 to ^^scertain the 
impact of horizontal resolution. 

Charts VIII-63 through VIII-66 contain the forecast 
error /change patterns at both sea level and 500 MBS. In 
addition, the actual change patterns are shown in Charts 
VIII-39 and VIII-60. T \ble VII-11 contains the actual change 
statistics. 

Table VII-12 contains model statistics for Run 
P22. At sea level, the 4.78-MB EtlS forecast change compares 
well to the 4.81-MB RM'S actual change, and the 3.61-MB RMS 
forecast error is much less than the actual change. At 500 
MBS, the 53.7-meter RMS forecast change compares well to the 
54.7-meter actual change, while the 36. 3-meter RMS error is 
much smaller than the actual change. 


2 . 


Run F18 (PECHCV) 


Charts through VIII-’136 contain the 

24-'hour sevi-level pressure, 500-MB height and 500-MB tempera- 
ture forecasts, respectively, using the 63 x 63 x 5 model 
version. Charts VIII-2 and VIII-6 may be used to verify the 
PS and 2500 forecasts, respectively. Charts VIII-137 through 
VIII-139 contain the 48-hour PS, 2500 and T500 forecasts, 
respectively. Verify the PS and 2500 forecasts using Charts 
VIII-3 and VIII-7, respectively. 

Cliarts VIII-67 through VIIJ-74 contain 'i.he forecast 
change/error patterns. Twenty- four hour changes/errors at 
both sea level and 500 i4BS are shown in Charts VIII-67 through 
VIII-7 0. Forty-eight hour changes/error’s are shown in Charts 
VIII-71 through VIII-74. These patterns may be compared to 
the actual changes shown in Charts VIII-59 through VIIl-62. 

Table VIl-13 contains model statistics for Run P18. 
At sea level, the 5.74-f'IB RMS forecast change is slightly 
less than the 6.01-MB actual change after 48 hours. Also, the 
4.16-MB RMS error shows good skill compared to the actual 
change. At 500 MBS, the statistics are consistent. The 
64.1-raeter RMS forecast change is slightly less than the 
67.7-meter RMS actual change. And, the 51.2-meter RMS error 
shows skill compared to the actual change . 
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3 . Run T4 (PECHFV) 

Charts VIII-140 through VIII-142 contain the 
24~hour PS,. Z500 and T500 forecasts, respectively, using 
the 63 X 63 X 10 model version. Charts VIII-2 and VIII-6 
may be used to verify the PS and Z5Q0 forecasts, respect- 
ively. Charts VlII-143 through VIII-145 contain the corres- 
ponding 48-hour forecasts. Verify the PS and Z500 forecasts 
using Charts VIII-3 and VI II -7, respectively. 

Charts VIII-75 through VIII-S2 contain the forecast 
change/error patterns. Tvjenty-four hour changes/errors at 
both sea level and 500 MBS are shown in Charts VIII-75 through 
VIII-78. Forty-eight hour changes/errors are shown in Charts 
VIII-79 through VIII-82. These patterns may be conipared to 
the actual changes shown in Charts VIII-59 through VIII-62. 

Table VIII-14 contains model statistics for Run 
T4. at sea level, the 5.81-MB RMS forecast change compares 
well to the 6.01-MB RMS actual change, while tlie 4.15-MB RMS 
error is much smaller than tlie actual change. At 500 MBS, 
the 64. 9 -meter RMS forecast change compares well to tlie 
67,7-meter RMS actual change, while the 50. 9-meter RMS 
forecast error shows good skill after 48 hours. 
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4. Run F20 (PECHCV) 

Charts VTII-146 through VIII-148 contain the 
24-hour PS, Z500 and T5Q0 forecasts, respectively, using 
the 63 X 63 X 5 model version and Production Series C 
analyses. Charts VIII-2 and VIII-6 may be used to verify 
the PS and Z500 forecasts, respectively. The forecasts ma.y 
be compared also to Charts VIII-134 through VIII-136 to 
ascertain the differences due to variations in the initial 
analyses. 

Charts VIII-149 through VIII-151 contain the 
48-hour PS, Z500 and T500 forecasts, respectively, for Run 
F20. Verify these forecasts against Charts VIII-3 and VIII-7, 
respectively. Compare these forecasts against VIII-137 
through VIII-139 to ascertain the impact of initial analyses 
on 48-hour forecasts. 

Charts VIlI-83 through VIII-90 contain the F20 
forecast change/error patterns. Compare these charts to the 
actual change patterns shown in Charts VXiI-59 through 
VIII-62. ■ ■ 

Table VIII-15 contains model statistics for Run 
P20. ht sea level,, the 4.26-tiB RJIS error shows skill compared 
to the 6,01-MB RllS actual change. At 500 MBS, the 56. 9-meter 
RMS error shows skill compared to the 67, 7-meter RtlS actual 
change. 
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5. 


Run (PECHFV) 


Charts VIlI-152 through VIII-154 contain tlie 
24~hour PS^ Z500 and T500 forecastSi- respectively, using 
the 63 X 63 X 10 model version Eind Production Series C 
analyses. Charts VI I I -2 and VIII-6 may be used to verify 
the PS and 2500 forecasts, respectively. The T6 forecasts 
may be compared also to Run F20 forecasts to ascertain 
vertical resolution effects, or to Run T4 forecasts to 
ascertain initial analysis effects. 

Charts VIlI-155 through VIII-157 contain the 
48-hour PS, S50Q and T500 forecasts, respectively, for Run 
T6. Verify these against Charts VIII-3 and VI I 1-7. Compari- 
sons to appropriate charts for Runs P20 and T4 may be made 
also. 

Charts VIXI-9r through VIII-9 8 contain the T6 
forecast change/error patterns- Compare these patterns to 
the actual change patterns shown in Charts VlII-59 through 
VIll-62. ; : 

Table VII-16 contains model statistics for Run T6. 

At sea level, the 4.13-MB RMS forecast erroi* shows skill com- 
pared to the 6.01 -mb RMS actual change. At 500 MBS, the 56.6- 
metef RMS forecast error shows skill compared to the 67. 7-meter 
RMS actual change. In both cases, the forecast change was only 
slightly less than tile actual change. 
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B. Scenario B: Foracas-bs/ Forecasl: Changes, Forecast 

Errors > ; 

. . . I ■ 

i ■ i 

■ ■ ■ :..J ■ 

All of the Scenario B forecasns were made from the same ^ . 

starting analyses. The 187 x 187 x 5 model version was 
initialized with values that were obtained by interpolation ; ; 

of 63 X 63 analysis model outputs. 

1. Run F24 (PEFHCV) 

f i 

Charts VIII-158 through VIII-160 contain the i ! 

24-hour PS;. Z500 and T500 forecasts, respectively, using | V 

the 187 X 187 x 5 model version and interpolated Production 
Series C 63 x 63 analyses. Charts VIII-14 and VII 1-18 may 
be used to verify the PS and 3500 forecasts, respectively. 

These forecasts may be compared also to Run P13 forecasts, 
shown as Charts VIII-161 through VIII-163, to assess the 
effects of horizontal resolution. 

Charts VIII-103 through VIII-106 contain the ,i 

forecast changs/error patterns at both sea level and 500 MBS. 

In addition, the actual change patterns are shown as Charts 
VlII-99 through VIII-1 02. Table VII-17 contains the actual 
change statistics. 

Table VII-18 contains the model statistics for 

... fci. ^ 

Run F24. At sea level, both the forecast change and the 

. . 

forecast error exceed the actual change. The RMS forecast 
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error was 3.89 MBS while the RMS actual change was 3.34 MBS. 
At 500 MBS, the 38.2-meter RMS error is only slightly 
smaller than the 38.6-meter R24S actual change. 

2. Run FI 9 (PECHCV) 

Charts VI 11-161 through VIII-163 contain the 
24-hour PS, Z500 and T500 forecasts, respectively, using the 
63 X 63 X 5 model version. Verify the PS and Z500 forecasts 
against Charts VIII-14 and VIII-18, respectively. Compare 
the F19 forecasts to the T5 forecasts to assess the effect 
of verticaJ. resolution. Charts VIll-164 through VIII-166 
contain the corresponding 48-hour forecasts for Run P19. 
Verify these against Charts VIII-15 and VIII-19, as 
appropriate . 

Charts VIlI-107 through VIII-114 contain the 
24-hour and 48-hour forecast change/error patterns. These 
patterns may be compared to the actual changes shown in 
Charts VIII-99 through VIII-102. Table VII-17 contains the 
actual change statistics. 

Table VII-19 contains the model statistics for 
Run F19, At both sea level and 50C MBS, F19 shows slight 
skill after 48 hours. At sea level, the 3.83-MB RMS forecast 
error is slightly smaller than the 4. 3 5 -MB RMS actual change. 


At 500 MT;S, the 50.8-meter RjyiS forecast error is slightly 
smaller than the 52.3-meter RMS actual change. This is 
true in spite of the fact that the forecast error exceeded, 
the actual change during the first 24 hours of the forecast. 

..i 

3. Rtm T5 (PECHFV) 

Charts VII 1-167 through VIII-169 contain the 
24-hour PS, Z500 and T500 forecasts, respectively, using 
the 63 X 63 x 10 model version. Verify the PS and 2500 
forecasts against Charts VIII-14 and VIII-18, respectively. 
Charts VIII-170 through VIII-172 contain the corresponding 
48-hour PS, 2500 and T500 forecasts. Verify these against 
Charts VIII-15 and VIII-19, as appropriate. 

Charts VIII-115 through VIII-122 contain the 
24-hour and 48-hour forecast change/error patterns. These 
may be compared to the actual changes shown in Charts VI II- 9 9 
through VIII-102. Table VII-17 contains the actual change 
statistics . 

Table VII-20 contains the model statistics for 
Run T5. At sea level, the 3.34-MB RMS forecast error after 
24 hours exactly equals the actual change, but the 3,70-MB 
RMS error after 48 hours is smaller than the 4.35-MB RMS 
actual change. At 500 MBS, the forecast error is slightly 
smaller than the actual change. After 24 hours, the 38.0-meter 
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RMS error may be compared to the 38,6-meter RMS change. 

After 48 hours, the 50,4-meter RMS error may be compared to 
the 52, 3 -meter RMS change. 

In general, the models did not exhibit as much 
skill in Scenario B (the May case) as they did in Scenario A 
(the April case) . 

C- Scenario A; Comparative Perf 02 nnaiice Data 

1. Intensity of Pressure Systems at Sea Level 

Table VTI-l contains the 24-hour forecast central 
pressures (in millibars) using three models of differing 
resolution. Eight significant low-pressure systems are 
considered. The starting and verification pressures for 
each system are provided. The average algebraic error is 
also shown for each model version. 

Run P22, the 187 x 187 x 5 model, exhibits the best 
overall performance in terms of algebraic error (-0,8 MBS per 
system) , due in part to large compensating errors (+11 MBS in 
the Alaskan Gulf low and -10 MBS in the China low) . Both FI 8 
and T4 exhibited the same +6 MB average error, more or less. 
Indeed, the effects of doubling the Vertical resolution appear 
to be minimal. Both coarse-mesh models did not deepen several 
of the lows enough. 


V-ll 



2 . 


Performance Statistics 


Tables VII-3 and VII-4 contain comparative per- 
formance statistics at 500 MBS and sea level/ respectively. 

In Table VII-3 , Scenario A, Run F22 has the best 
24-hour performance, with an RMS height error of only 36.3 
meters. P22 also has an RMS forecast change of 53.7 meters, 
vzhich compares very favorably to an actual PJys change of 
54.7 meters. At 48 hours, the 63 x 63 models (F18 and T4) 
exhibit RMS errors which are quite comparable to one another - 

In Table VII-4, Scenario A, Run F22 has a slightly 
larger RMS error (3,61 fIBS) than the 63 x 63 models. It did, 
however, produce an RMS forecast change which compares very 
well to the RMS actual change. Both 63 x 63 models tended to 
produce small error as a consequence of a less-courageous 
(smoother) forecast. 

3, Kinetic Energy (K E) 

Table VII-5 contains the layer-mean kinetic energy 
24-hour forecast changes (expressed in percent) for models of 
varying resolution. Consider the Scenario A part. Run F22 
produced large increases in KE in the lower (moisture bearing) 
levels, as compared to the 63 x 63 models which both predicted 
KE losses at all levels. Figure VI-1 shows the percent fore- 
cast change in KE by model level. 
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Figures VI-2 through VI-4 contain KE 24-hour fore- 
cast cross-sections (latitude vs. model level ) , in energy 
units. Figure VI-2 pertains to Run T4, which exhibits a 
maximum of 18 units hear 300 MBS at 40®N, Figure VI-3 per- 
tains to Run P18, and exhibits a 16+ unit maximum near 400 MBS 
at 40°N, Altitude differences are due to vertical resolution. 
Figure VI-4 exhibits a stronger maximum (22+ units) near 
300 MBS and 40°H. Thus^ the 187 x 187 x 5 model and 
63 X 63 X 5 model exhibit comparable KE core locations , but 
the 187 X 187 model produces a much greater core magnitude. 

Figures VI-5 through VI-9 show the time variations 
of KE by model/model level for Scenario A. In Figure VI-S, 

Run F22 exhibits a continuing growth with time at sigma=0.9 
level, as compared to Runs PlS and T4 which decrease after 
reaching a 6-hour- .naximum. In Figure VI-6, one sees -fhe same 
type of behavior at sigma=0.7 level. In Figure VI-7, Run F22 
exhibits some oscillatory behavior without the large net 
increase observed at lower levels. Runs F18 and T4 show net 
decreases. In Figure VI-8, one can see that all three models 
behave in like manner — wi-th a large initial decrease 
followed by oscillatory beihavior. (This appears to. be 
typical at the sigraa=0,3 level. ) Finally, the somewhat less- 
remarkable behavior at the sigma=0.1 level is shown in Figure 
VI- 9. 
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4. 


Square Vorticity 


Figiire VI-10 shows the square vorticity 24-hour 
forecast changes (in percent) by model version and level. 
PEFHCV (Run F22) produced large increases in the lower 
levels compared to PECHCV (F18) and PECHFV (T4) . (Table 
VII-6 contains the source percentages.) 

Figures VI-11 through VI-15 contain the time 
variations of square vorticity at each level for Scenario A. 
Run F22 (187 x 187 model) generates much more vorticity than 
the tv70 63 X 63 models in the lowest three levels. In Figure 
VI-14, one can observe a decrease in square vorticity in the 
first six hours (in all models) at the sigma=0.3 level. Com- 
pare this to the KE behavior at this level shown previously 
in Figure VI-8. 


5. Square Divergence 

Figures VI-16 through VI-20 contain the time vari- 
ations of square divergence in Scenario A at the five model 
levels, respectively. Rim F22 (the 187 x 187 model) behaves 
quite differently from 63 x 63 models, especially at low 
levels. The effects of vertical resolution am? horizontal 
resolution appear to be (more) alike at higher levels. 
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6 . 


Temperattire 


Figures VI -41 through VI-43 contain 24-hour fore- 
cast temperature-change (°C) cross-sections (latitude versus 
model level) for Scenario A. forecasts. Figure VI-41, which 
pertains to Run T4, and Figure VI-42, which pertains to Run 
F18, are in qualitative agreement. Figure VI-43, which 
pertains to Run F22, shows a greater (positive) forecast 
change (near 500 MBS at 40°N) than the 63 x 63 model versions. 

7. Precipitation 

Figure VI-47 shows the average precipitation 
(centimeters) per grid point for three model versions in 
Scenario A. The source numbers are contained in Table VII-7 . 
The effects of vertical resolution changes are small compared 
to the effects of horizontal resolution changes. 

D. Scenario B; Comparative Performance Data 

1, Intensity of Pressure Systems at Sea Leve l 

Table VII-2 contains the twenty-four hour forecast 
central pressures (in millibars) using three models of dif- 
fering resolution. Six significant low-pressure centers are 
considered in this scenario. The starting and verification 
pressures for each system are provided. The average algebraic 
error for each model version is also tabulated. 
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Runs T5 and F19, the 63 x 63 models, exhibit the 
best overall performance. As in Table VII~1, the effects of 
vertical resolution appear to be minimal. In contrast. Run 
P24 (the 187 x 187 x 5 model) produced lows which were 
slightly too deep ( 'vS MBS on the average) . 

2. Performance Statistics 

Table VII~3 and VII-4 contain comparative per- 
fomance statistics at 300 DIBS and sea level, respectively. 

In Table VII-3, Scenario B, all three model 
versions perform in the same manner. The R14S error at 
500 MBS is equal to the RMS change (no skill compared to 
persistence) at 24 hours, but slightly less at 48 hours. 

In Table VII-4, Scenario B, T5 exhibits the smallest 
RMS error, but mainly as a consequence of a less-courageous 
(smoother) forecast. Run P24 does not exhibit skill compared 
to persistence at 24 hours. It predicted too much change. 

The 63 X 63 models show slight skill at 48 hours, however. 

3 . Kinetic Energy (KE) 

Table VI I -5 contains the layer-mean kinetic energy 
24-hour forecast changes (expressed in percent) . Consider 
the part for Scenario B, F24 produced a large increase in 
KE (27.8%) at sigma=0.9 level, compared to KE decreases in 


the 63 X 63 models. At other levels, F24 produced smaller 
decreases in KE than the 63 x 63 models. Figure VI-21 
shows the percent forecast change in KE by model level. 

Figures VI-22 through VI-24 contain KE 24-hour 
forecast cross-sections {latitude versus model level) in 
energy units. Figure VI-22 shows the Run T5 cross-section 
with a core magnitude of 12+ tmits near 300 MBS at 40°N. 
Figure VI-23 shows the Run P19 cross-section with a core 
magnitude of 12+ units near 350 MBS and 45°N. Figure VI-24 
shows the Run F24 section with a core magnitude of 12+ units 
near 350 MBS at 45°N. In this May case, the core magnitudes 
are quite comparable. The ten-level model places the core 
slightly higher in altitude than the five-level models. 

Figures VI-25 through VI-29 contain the time vari- 
ations of kinetic energy by model/model level for Scenario B. 
In many ways, the KE behavior is comparable to the behavior 
described for Scenario A earlier. In Figure VI-25, one 
observes that the sigma=0.9 level net increase by Run P24 
(187 X 187 model) is much xess than that shown in Figure 
VI-5 (by Run F22) . The modal value at forecast hour six 
remains for the 63 x 63 models. In Figure VI-26, however, 
the F24 behavior at s5gma=0.7 level is much better than F22 
behavior shown previously in Figure VI-6. In Figure VI-28, 
one sees the (same) large decreases in the first six hours at 
the sigma==0.3 level. 
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4. 


Square Vorticity 


Figure VI-30 shows the square vorticity 24-hour 
forecast changes (in percent ) , by model version and level. 
PEFHCV (Run F24) produced a .large increase at the sigma=0.9 
level. Compare this to the behavior of the 63 x 63 models 
shown. (Table VII-6 contains the source percentages.) 

Figures VI-31 through VI-35 contain the time 
variations of square vorticity at each level for Scenario 
B. Run F24 (the 187 x 187 model) behavior is quite different 
from the 63 x 63 models. 

5. Square Divergence 

Figures VI-36 through VI -40 contain the time 
variations of square divergence in Scenario B at the five 
model levels, respectively. R\jn F24 behaves quite differ- 
ently from 63 x 63 models, especially at low levels. The 
effects of vertical resolution and horizontal resolution 
appear to be comparable at higher levels. 

6. Temperature 

Figures VI -44 through VI-46 contain 24-hour fore- 
cast temperature-change (“C) cross-sections (latitude versus 
model level) for Scenario B forecasts. Run T5 in Figure VI-44 
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shows qualitative agreement with Run F19 in Figure VI-45, 
but it is not as smooth. Run F24 in Figure VI-46 corresponds 
very well with F18 but tends to show larger amplitudes near 
35“N. 

7. Precipitation 

Figure VI-48 shows the average precipitation 
(centimeters) per grid point for three model versions in 
Scenario B. The source numbers are shown in Table VII-7. 

As in Scenario a!, the effects of vertical resolution changes 
are demonstrated to be small compared to horizontal reso- 
lution changes. 

E. Effects of Changes in Model Vertical Resolution 

on Sea- Level Pressure and 500-I4B Height Forecasts 

1. Scenario A Comparisons 

Three comparisons can be made (in Scenario A) 
between five- layer and ten-layer models: (1) Runs F19 and T4 
which were both initialized with Test Series B analyses; (2) 
Runs F20 and T6 which were both initialized with Production 
Series C analyses; and (3) , an additional set of forecasts 
which were initialized with analyses using FNWC first-guess 
fields. (The comparison between Run F22 and its ten-layer 
counterpart will be covered in a separate volume.) 
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Chaxts VIII-37 and VIII‘-38 show the differences 
between Runs P20 (63 x 63 x 5) and T6 (63 x 63 x 10) at sea 
level and 500 MBS, respectively, (Chart VIII-37 has a 2 -MB 
contour interval; and Chart' VIII-38 has a 30-meter contour 
interval.) Table VII-21 contains statistical measures of the 
differences between P20 and T6 at the 24-hour and 48-hour 
forecast times at both sea level and 500 MBS. The RMS 
pressure difference is 0.96 MBS at hours 24 and 48 (no 
additional differences after the first 24 hours). At 500 MBS, 
the RMS height difference increases from 8-34 meters at hour 
24 to 9.32 meters at hour 48. Prom Chart VIII-37, note that 
a maximum 24-hour difference of 5.15 t4BS is located 35N/160E. 

Charts VIII-38 and Vj'II-39 show the differences 
between Runs P18 and T4. Table VII-22 contains statistical 
measures of the differences. In this case, the RMS pressure 
difference is 0.92 MBS at hour 24 and 1,02 MBS at hour 48. 

The maximum difference, however, is 7.29 MBS near 35N/160E. 

At 500 MBS, the RMS height difference increases from 8.19 
meters at hour 24 to 9.79 meters at hour 48. 

Charts VIII-41 and VIII-42 show the differences 
between a five-layer and ten-layer forecast, both of which 
were initialized with analyses that started with PNWC first- 
guess fields. Prom visual appearances, the differences appear 
to be slightly smaller than those observed in the previous cases. 
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These differences appear to be small in both 
magnitude and areal extent. They can be compared to differ- 
ences between a 187 x 187 forecast and a 63 x 63 forecast 
which were approximately twice as large in magnitude. 
Finally, one can also assess the forecast differences 
brought about by differences in the initialization. In an 
RMS sense, such differences appear to be small compared to 
the differences brought about by changing the vertical reso- 
lution. 


2. Scenario B Comparison 

Charts VIII-43 and VIII-44 show the differences 
between Runs F19 and T5. Table VII-23 contains statistical 
measures of the differences at sea level and 500 MBS. The 
RMS pressure difference of 0.76 MBS at hour 24 increases to 
0.82 MBS at hour 48. The maximum 24-hour difference is 
+4.88 MBS south of Japan. At 500 I4BS, the RMS height differ- 
ence of 7.27 meters at hour 24 increases to 8.66 meters at 
hour 48. NOTE: The effects of differing vertical resolution 
on model energetics, precipita;-ion and forecast accuracy 
were discussed earlier. Suffice to say that the differences 
due to vertical resolution were consistently small {or even 
negligible) compared to the differences due to horizontal 
resolution. Refer to Table VII-7, for example, which shows 
the precipitation amounts from each model. 
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F- Effects of Changes in Model Horizontal Resolution 

on Sea-Level Pressure and 5Q0-MB Height Forecasts 

1. Scenario A Comparison 

Charts VIlI-45 and VIII-46 show the differences 
between Runs P22 (187 x 187 x 5) and F18 (63 x 63 x 5) at 
sea level and 500 liBS, respectively. (The sea-level pressure 
contour interval is 2 MBS. The 500-MB contour interval is 
30 meters.) Table VII-9 contains statistical measures of 
the differences. The RMS pressure difference is 2.01 MBS, 
v;ith a raaximina difference of -16.5 MBS near Japan. At 500 
I4BS, the RMS height difference is 14.1 meters, with a maximum 
difference of -148.0 meters near Japan, These differences 
are "of meteorological scale”. They tend to contour nicely. 
The differences in the subtropics are small. 

2, Scenario B Comparison 

Charts VIII-47 and VIII-48 show the differences 
between Runs F24 (187 x 187 x 5) and P19 (63 x 63 x 5) at 
sea level and 500 MBS, respectively. Table VII-10 contains 
statistical measures of the differences. The RMS pressure 
difference is 1,58 MBS, with a maximum difference of -9.6 MBS, 
At 500 MBS, the ]^IS height difference is 9.8 meters, with a 
maximum difference of -69,6 meters. 
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NOTE: Differences in model energetics, precipitation, and 
forecast accuracy as a consequence of changes in horizontal 
resolution were discussed earlier. Such differences are large 
compared to other causes of differences, generally. 
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FIGURES 
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FIGURE VI- 1; KINETIC ENERGY 24-HOUR FORECAST 

CHANGE (%) , BY MODEL VERSION 
AND LEVEL. SCENARIO A. 


PAQg 


m 








psmJNAL PAGE IS 
pS POOR QUALITY 























IT 


I 


U 


i 

1 


I 



FIGURE. VI- 7 f 


of: kiiietic : energy: for 
:;^hbeN;:modeN::^ ivc' siGMa=o.'5-^ 

SGEN^ 






KINETIC ENERGY (xlO) 


SIGMA=0.3 



16 12 18 24 


TIME VARIATION OF KINETIC ENERGY FOR 
THREE MODEL VERSIONS AT SIGI4A=0 . 3 
LEVEL. SCENARIO A. 
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FIGURE VI-13: TIME VAHIATION OF SQUARE VORTICITY FOR THREE 

MODEL VERSIONS AT SIGMA=0.5 LEVEL. 

SCENARIO A. 
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FIGURE VI- 15: TIME VARIATION OF SQUARE VORTICITY FOR THREE 

MODEL VERSIONS AT SIGMA=0.1 LEVEL. 

SCENARIO A. 
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FI-GURE VI- 17: TIME VARIATION OF SQUARE DIVERGENCE FOR THREE 

MODEL VERSIONS AT SIGMA=0.7 LEVEL. 

SCENARIO A. 
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FIGURE VI- 18: TIME VARIATION OF SQUARE DIVERGENCE FOR THREE 

MODEL VERSIONS AT SIGMA=0.5 LEVEL. 

SCENARIO A. 
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FIGURE VI-19; TIME VARIATION OF SQUARE DIVERGENCE FOR THREE 

MODEL VERSIONS AT SIGMA=0.3 LEVEL. 

SCENARIO A. 
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FIGURE VI- 20: TIME VARIATION OF SQUARE DIVERGENCE FOR THREE 

MODEL VERSIONS AT SIGMA=0.1 LEVEL. 

SCENARIO A. 
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h FIGURE VI-21: KINETIC ENERGY 24-HOUR FORECAST 

CHANGE (%) , BY MODEL VERSION 

I AND LEVEL. SCEN2LRIO B. 
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VI-23: KINETIC ENERGY (xlO) 24 -HOUR FORECAST CROSS-SECTION 

RUN F19. MODEL PECHCV. SCENARIO B. 
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FIGUEE VI-24: KINETIC ENERGY (xlO) 24-HOUR FORECAST CROSS-SECTION. 

RUN F24. MODEL PEPHCV. SCENARIO B. 
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VI-30: SQUARE VORTICITY 24-HOUR FORECAST 

CHANGE (%) , BY MODEL VERSION AND 
LEVEL. SCENARIO B. 
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FIGURE VI- 31 


TIME VARIATION OF SQUARE VORTICITY FOR THREE 
MODEL VERSIONS AT SIGMA=0.9 LEVEL. 

SCENARIO B. 
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FIGURE VI-40: TIME VARIATION OP SQUARE DIVERGENCE FOR THREE 

MODEL VERSIONS AT SIGMA=0*1 LEVEL. 

SCENARIO B. 





L&TITUDE 


FIGURE 


VI-41: 


TEMPERATURE 24-HOUR FORECAST CHANGE («C) , 
BY LATITUDE AND LEVEL. RUN T4. MODEL 
PECHPV. SCENARIO A. 


»«3E is 

















LATITUDE 


EIGDRE 


VI-45; TEMPERATURE 24-HOUR FORECAST CHANGE (“C) , BY LATITUDE 
AND LEVEL. RUN F19. MODEL PECHCV. SCENARIO B. 







LATITUDE 


nij:: 


FIGURE VI-46: TEMPERATURE 24-HOUR FORECAST CHANGE (®C) , BY LATITUDE AND 

LEVEL. RUN F24 . MODEL PEPHCV. SCENARIO B. 
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TABLE VII- 1 


24-HOUR FORECAST CENTRAL PRESSURES* (MBS) USING 
THREE MODELS OF DIFFERING RESOLUTION. 

SCENARIO A. 




— ^ ^ ^ ■ q 

LOW-PRESSURE CENTERS 

AVERAGI 

algebra: 

ERROR 

(MBS) 

SGEN^IO 

RUN 

CHINA 

EAST OF 
KAMCHATKA' 

ALASKAN 

GULF 

WEST 

ATLANTIC 

MED 

TEXAS - 
COLORADO 

JAMES 

BAY 

NOVAYA 

ZEMYLA 


F22: 

975 

1001 

991 

983 

■ 997 

' 

999 

994 

997 

-0,8 

' A 

P18 

986 

1009 

996 

990 

1001 

1003 

1002 

1000 

5.5 


T4 

986 

1009 

996 

990 

1002 

1004 

1004 

1000 

6.0 

STARTING 

ANALYSIS 

1002 

1003 

1003 

988 

1 

1004 

1005 

998 

992 


UERIFICATION 

ANALYSIS 

985 

998 

980 

986 

999 

1003 

1002 

990 V 



* FROM 1200 Z, 22 APRIL 1976 


STARTING 

ANALYSIS 


VERIFICATION 

ANALYSIS 


* FROM 1200Z, 20 MAY 197S 




























































TABLE VI I- 4 


PERPOEMANCE STATISTICS FOR FORECAST MODELS OF 
DIFFERING RESOLUTION: SEA-LEVEL PRESSURE (MBS) 


RMS DIFFERENCES 


2 4 -HOUR 


4 8 -HOUR 


ECENARIG 


RUN 


MODEL 


ACTUAL 

CHANGE 


FORECAST 

CHANGE 


ERROR 


ACTUAL 

CHANGE 


FORECAST 

CHANGE 


ERT 


A' 


P18 

T4 

F22 


PECHCV 

PECHFV 

PEPHCV 


4.81 


4.19 

4.21 

4.78 


3.54 

3.47 

3.61 


6.01 


5,74 

5.81 


4. 

4. 


B 


FI 9 

T5 

F24 


PECHCV 

PECHFV 

PEFHCV 


3.34 


3.14 

2.99 

3.66 


3.55 

3.34 

3.89 


4.35 


4.27 

4.18 


3.83 

3.70 



























TABLE VII-Ss KINETIC ENERGY 24-HOUR FORECAST CHANGES 
FOR MODELS OF VARYING RESOLUTION AND 
TWO SCENARIOS (GIVEN IN PERCENT) . 


SIGMA LEVEL 


SCENARIO 

MODEL 

0.9 

n , . 1. 1 _ 

0.7 

: 1 

0.5 

0.3 

0.1 


PECHCV 

(F18) 

“4.9 

“13.1 

-10.9 

-19.5 

-16.9 

A '■ 

PECHFV 

(T4) 

-8.4 

-12.5 

-13 . 8 

-18.4 

-11.2 


PEFHCV 

(F22) 

52.8 

24.3 

6.4 

-9.8 

-13.7 


PECHCV 

(P19) 

-13.7 

-27.6 

-23.7 

-23.5 

-29.4 

B 

PECHFV 

{T5) 

“22.3 

-30.6 

-26.8 

-24.9 

-19.3 


PEFHCV 

(F24) 

27.8 

-3.2 

-6.6 

“15.2 

-22.2 






















TABLE VII- 6: SQUARE VORTICITY 24-HOUR FORECAST CHANGES, 

FOR MODELS OF VARYING RESOLUTION AND TWO 
SCENARIOS (GIVEN IN PERCENT) 


SIGMA LEVEL 


SCENARIO 

MODEL 


0.9 

0.7 

0.5 

0.3 

0.1 


PECHCV 

(F18) 


-12.2 

-43.2 

-37.1 

-37.4 

-43.4 

A 

PECHFV 

(T4) 


-25.3 

-44.2 

-44.3 

-33.5 

-21.0 


PEFHCV 

(F22) 


127.7 

28.5 

12.8 

-5.9 

-28.0 


PECHCV 

(F19) 


-13.9 

-48.1 ! 

-46.1 

-36.6 

-43.1 

B 

PECHFV 

(T5) 


-31.6 

-52.3 

-50.3 

-40.1 

-28.0 


PEFHCV 

(F24) 


115.7 

12.8 

-1.6 

-13.9 

-25.6 



V y 

































TABLE VII- 8 


STATISTICAL DIPPEKENCES BETVTEEN 
187 X 187 AND 63 x 63 ANALYSES. 


A. SEA-LEVEL PRESSURE (MBS) 


SCENARIO 

DIFFERENCE STATISTICS 

RMS 

MEAN 

MAX 

MIN 

a1 

0.35 

-Q.03 

2.27 

-3.59 

b2 

■0.37 

1 

0 

1 

o 

2.36 

-3.15 


B. 500 MB HEIGHT (METERS) 


SCENARIO 

DIFFERENCE STATISTICS 

RMS 

MEAN 

MAX 

MIN 

A 

7.46 

-1.67 

31.3 

-45.3 

B 

9.03 

-1.84 

33.2 

-98,8 


1. A = 1200Z, 22 April 1976 

2. B = 1200Z, 20 May 1976 

































TABLE VII- 9 


STATISTICAL DIFPEEENCES BETWEEN 
187 X 187 FORECAST AND 63 x 63 
FORECAST. SGENJ^IO A. 


A. SEA-LEVEL PRESSURE (MBS) 


HOUR 

DIFFERENCE STATISTICS 

RMS 

MEAN 

MAX 

MIN 

24 

2.01 

-0.81 

5.71 

-16.5 


B. 500 MB HEIGHT (METERS) 


HOUR 

DIFFERENCE STATISTICS 

RMS 

MEAN 

MAX 

KIN 

24 

14.1 

0.20 

89.3 

-148.0 























: TABLE ^X-'IO: iaTATISTIGALrBIEFEEffiNeES' B 

' -r;: ,:eoisgast 63 x 63. , 


^ F:^SSDRE (MBS) 


HOUR r 


Bias' 

•MEAN 

MAX 

MIN : 

24 ■■ 

iVsa 

;. -0.68 

4.50 

"y . ■: ; 

-9.6 



























TABLE VII-11: ACTUAL CHANGES (PINAL MINUS 

INITIAL ANALYSIS) . 

SCENARIO A. 


A. SEA-LEVEL PRESSURE (MBS) 


HOUR 


DIPPERENCE 

STATISTICS 


RMS 

. MEAN 

MAX 

MIN 


to 

H 

4,81 

0.13 


21.0 

-27.4 

48^ 

6.01 

-0.07 


25.3 

-34.6 


B. 500 MB HEIGHT (METERS) 


rjnTT'D 


DIPPERENCE 

STATISTICS 


IIVJUXV 







RMS 

MEAN 


MAX 

MIN 

24 

54.7 

3.1 


242.0 

-397.0 

48 

67.7 

-2.59 


359.0 

-296.0 


1. 4-23-76 MINUS 4-22-76 

2. 4-24-76 MINUS 4-22-76 
















TABLE VII- 12: MODEL STATISTICS. RON MODEL 

PEEHGV. SCENARIO A. TEST SERIES 
B ANALYSES. 


A. SEA-LEVEL PRESSURE CMBS) 


HOUR 

TYPE 


.* ■; ' , DIEEERENGE ^ 

STATISTICS 

— 


MAX 

MIN 


FORECAST CHANGE 


4; 78 

0.06 ■ 

16.9 

^38.0 

24 

ERROR 


3'. 61 

-0.42 

20.4 

-20.7 


ACTUAL CHANGE • 


4.81 

'0.13 

21.0 

-27.4 • 


B, 500 MB HEIGHT (METERS) 




























TABLE VII-13: MODEL STATISTICS. RUN Fl8 , MODEL 

:?ECHCV. SCENARIO A. TEST SERIES 
B ANALYSES. 


A. SEA-LEVEL PRESSURE (MBS) 


HOUR 

TYPE 

DIFFERENCE STATISTICS 

RMS 

MEAN 

MAX 

MIN 

24 

FORECAST CHANGE 
ERROR 

4.19 

3.54 

m 

18.2 

22,4 

-24.8 

-13.8 

48 

FORECAST CHANGE 
ERROR 

ACTUAL CHANGE 


■ 

22.2 

19.6 

25.3 

-37.2 

-20.1 

-34.6 

B. 500 

MB HEIGHT (METERS) 

HOUR 

TYPE 

DIFFERENCE STATISTICS 

RMS 

MEAN 

MAX 

MIN 

24 

FORECAST CHANGE 
ERROR 

47.1 

39.1 

23.3 

6.3 

173.0 

241.0 

-298.0 

-189.0 

48 

FORECAST CHANGE 
ERROR 

ACTUAL CHANGE 

64.1 

,51.2 

67,7 

24.9 

13.5 

-2.59 

278.0 

294.0 

359.0 

-329.0 

-194.0 

-296.0 
































TABLE'- :• MODEL.vSTATISTICS ■..RUN. 

vV.-: ; : . \ i'k^scenariq . aw:--^; test series- 

• B ANALYSES.. 


,Aw:'v-" sea-level, pressure- .(mbs-) 



TYPE- - -. . 


DIFFERENCE STATiSTIGS 


•RMS’ ' -S-v mAN 


FpRECAST CHANGE 

4:.:2Iv>.: 

;>',|5RRpR . 

i-i ,3-'47w; 

FORECAST CHANGE 

■: PsAai. ■ 

ERROR 

' : 44 . 5 ,..- 

ACTUAL Cm^GE; 

■6..01. 



0.59 . 


v2i.P-;:' 



■ r0i 


''.pi2;^ 

9 

CO 

: 



, - 34 I 



B. .. ..500...MB HEIGHT. (METERS). . 



































TABLE VII~15: MODEL STATISTICS. RUN F20. MODEL 

PECHCV. SCENARIO A. PRODUCTION 
SERIES C ANALYSES. 


A. SEA-LEVEL PRESSURE (MBS) 


HOUR 

TYPE 

DIFFERENCE STATISTICS 

RMS 

MEAN 


MIN 

24 

FORECAST CHANGE 
EPvROR 

4.01 

3,84 

m 

14.1 

22.5 

-27.2 

-18.1 

48 

FORECAST CHANGE 
ERROR 

ACTUAL CHANGE 

5.63 

4.26 

6.01 

-0.43 

-0.71 

-0.07 

20.2 

20.2 

25.3 

-31.6 

-25.7 

-34.6 


B, 500 MB HEIGHT (METERS) 


HOUR 

TYPE 

DIFFERENCE STATISTICS 

RMS 

MEAN 

MAX 

MIN 

24 

FORECAST CHANGE 
ERROR 

40.5 

48.2 

-5.86 

-22.8 

148.0 

255.0 

-316.0 

-220.0 

48 

FORECAST CHANGE 
ERROR 

ACTUAL CHANGE 

60.4 

,56.9 

67.7 

-3.99 

-15.3 

-2.59 

228.0 

310.0 

359.0 

-314.0 

-233.0 

-296.0 






































TABLE VII~16: MODEL STATISTICS . RtM MODEL 

PEGHFV, SCENARIO A. PRODUCTION 
SERIES C ANALYSES. ““ 


A. SEA-LEVEL PRESSURE (MBS) 


HOUR 

TYPE 

DIFFERENCE 

STATISTICS 

RMS 

MEAN 

• 

MAX 

MIN 

LZ ^ . _ : . ' n 

24 

FORECAST CH21NGE 

3.94 

-0,07 

14.4 

-23.6 


ERROR 

3.62 

-0.54 

21.8 

-16.3 

48 

FORECAST CHANGE 

5 . 62 

-0.11. 

21 .a 

: -32.3 


ERROR 

4.13 

-.39 

19.3 

-19.8 


ACTUAL CHANGE 

6,01 

-0.07 

2:5.3 

-34.6 


B. 500 MB HEIGHT (METERS) 


HOUR 

TYPE 

DIFFERENCE 

STATISTICS 

RMS 

MEAN 

MAX 

MIN 

■■ - ■■ ^ 1 

■ 24 

FORECAST' CHANGE 



40.9 

-4.61 

150.0 

-318 „0 


ERROR 

46.7 

-21.6 

254.0 

; -214.0 

48 ■■ 

FORECAST CHANGE 

60.6 

-4.01 

291.0 

-296.0 


error; 

.56.6 

-15.4 

305.0 

^ -234.0 


ACTUAL CHANGE 

67.7 

r2.S9 

359.0 

; -296.:o 






TABLE VII-17: ACTUAL CHANGES (FINAL MINUS 

INITIAL ANALYSIS) . 

SCENARIO B. 


A. sea-levEl pressure (mbs) 



DIFFERENCE STATISTICS 

JU^ 

1 

Riys 

MEAN 

MAX 

MIN 

to 

3.?4 

-0.07 

17.8 

1 

H 

• 

48^ 

4.35 

-0.18 

26.0 

-25.3 


B. 500 MB HEIGHT (METERS) 


HOUR I- 

RMS 

24 38.6 


DIFFERENCE STATISTICS 

MEAN MAX MIN 

-1.51) 201.0 -224.0 


48 52.3 2,87 332.0 -338.0 


5-21-76 MINUS 5-20-76 
5-22-76 MINUS 5-20-76 













0?ABLE VII-18: MODEL STATISTICS. RUN F24 . MODEL 
PEPHCV. SCENARIO B. PRODUCTION 
SERIES C ANALYSES. 


A, SEA-LEVEL PRESSURE (MBS) 


HOUR 

1 

TYPE 

•; DIFFERENCE STATISTICS 

RMS 

MEAN 

MAX 

MIN 




3.66 

-0.86 

15.2 

-17.5 

3.89 

-1.06 

19.6 

-19.1 

3.34 

-0.07 

17.8 

-17.9 


24 


FORECAST CHANGE 
ERROR 

ACTUAL CHANGE 


B. 500 MB HEIGHT (METERS) 




























TABLE VIT-19: MODEL STATISTICS. RUN E19 . MODEL 

RECHCV. SCBI^ARIP B« PRODUCTION 
SERIES C ANALYSES . 


A, SEA-LEVEL PRESSURE (MBS) 


HOUR 


TYPE 


DIFFERENCE STliTX-STICS 


MEAN 


24 

FORECAST CHANGE 
ERROR 

CO 

FORECAST CHANGE 
ERROR 

ACTUAL CHANGE 


3.14 

3,55 


4.27 

3.83 

4.35 


■G.14 

-0.23 

-0.18 


13.7 

20.5 


23.1 
15.6 
26. C 


-18.9 

-14.6 

-25.3 


500 MB HEIGHT (METERS) 


HOUR 

TYPE 

DIFFERENCE 

STATISTICS 

RMS 

MEAN 

MAX 

MIN 



24 

FORECAST CHANGE 

35.4 

2.44 

156 .0 

-216.0 

ERROR 

39.3 

-13. 8 

183. 0 

-140.3 

: i 

48 

j 

FORECAST CHANGE 

52.2 

-0.66 

238.0 

-194.0 

ERROR 

.50.8 

-16.4 

256 . 0 

-240.0 


ACTUAL CHANGE 

52.3 

2,87 

332.0 

-338.0 






































TABLE VII-20: MODEL STATISTICS. RUN MODEL 

PECHFV . SCENARIO B. PRODUCTION 
SERIES C ANALYSES. 


A. SEA-LEVEL PRESSURE (MBS) 


HOUR 

TYPE 

DIFFERENCE 

STATISTICS 

RMS 

MEAN 

MAX 

MIN 



n 

24 

FORECAST CHANGE 

2.98 

1 

o 

■ 

o 

13.7 

-13.8 

ERROR 

3.34 

-0.24 

20.3 

-14.6 

48 

FORECAST CHANGE 

4.18 

-0.10 

24,0 

-17.1 

ERROR 

3,70 

-0.19 

14.7 

-14.4 


ACTUAL CHANGE 

4.35 

-0.18 

26.0 

-25.3 


B. 500 MB HEIGHT (METERS) 


HOUR 

TYPE 

DIFFERENCE 

STATISTICS 

RMS 

MEAN 

MAX 

MIN 

--- 1 

24 

FORECAST CHANGE 

35.0 

-3.56 

146.0 

-222.0 


ERROR 

38.0 

-14.90 

182.0 

-138.0 

48 

FORECAST CHANGE 

52.7 

-3.19 

248.0 

-197.0 


ERROR 

.50.4 

-18.90 

246.0 

-246.0 


ACTUAL CHANGE 

52.3 

2.87 

332.0 

-338.0 




































I-2X: PECHFV POHSa^ST MINUS PEGHCV 

FOEECAS'E* SCENARIO A. PRO- 
DtrCilQN SERIES C ANALYSES. 


A, SEA-LEVEL PRESSURE (MBS) 


DIEEERBNCE , STATISTICS 

RMS 

..MEAN 

MAY 

MIN 

0.S6 

0.40 

5 .15 

-4.24 

0.96 

0.32 

5.90 

-4.54 


300 MB HEIGHT (METERS) 


DIPEERENCE STATISTICS 


1,25 

- 0.02 


MAX 

MIN 

38.4 

67.0 

“42.4 

“38.0 
























TABLE VII-22: PEGHPV FORECAST MINUS PECHCV 

FORECAST, SCENARIO A. TEST 
SERIES B ANALYSES. 


A. SEA-LEVEL PRESSURE (MBS) 


HOUR 


DIFFERENCE 

STATISTICS 


RMS . 

MEAN 

MAX 

MIN 


24 


0.20 


7.29 

“3,77 

48 

1^1 

0 . 15 


9.35 

-8.03 


B. 500 MB HEIGHT (METERS) 


DIFFERENCE STATISTICS 


HOUR 


RMS 


MEAN 


MAX 


MIN 


24 


8.19 


0.64 


48 


S .79 


-0.52 



















TABLE VII- 23: PECHPV FORECAST MINDS PECHCV 

FORECAST. SCENARIO B. PRO- 
DUCTION SERIES C ANALYSES. 


A, SEA-LEVEL PRESSURE (MBS) 



DIFFERENCE STATISTICS 

rxuujtv 


^ 




RMS 

•MEAN 

MAX 

1 

MIN 

24 

0.76 

0.14 

CO 

CO 

• 

-2.50 

48 

0.82 

0.40 

4.97 

1 

• 

to 


B. SOO MB HEIGHT (METERS) 


HOUR 

DIFFERENCE STATISTICS 

RMS 

MEAN 

MAX 

MIN 


24 

7.27 

-1.12 

28.3 

-23.4 

48 

8.66 

-2.53 

51.5 

-42.2 
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